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DEVELOPMENT OF AN ANALYSIS FOR THE DETERMINATION OF 
COUPLED HELICOPTER ROTOR/CONTROL SYSTEM DYNAMIC RESPONSE 
PART I - ANALYSIS AND APPLICATIONS* 

By Lawrence R, Sutton and Stephen A. Rinehart 
Rochester Applied Science Associates, Inc. 


SUMMARY 


A theoretical analysis is developed for a coupled heli- 
copter rotor system to allow determination of the steady- 
state dynamic response behavior of helicopter rotor systems in 
steady-state forward flight or maneuvers. The effects of an 
anisotropically supported swashplate or gyroscope control system 
and a deformed free-wake on the rotor system dynamic response 
behavior are included in this analysis. 

The analysis involves the utilization of a combination of 
Laplace transform and transfer matrix techniques such that a set 
of Laplace transformed equations in matrix form result which in- 
clude interharmonic coupling of the blades due to the control 
system and blade aerodynamic coupling. The solution of the matrik 
equations results in the determination of the harmonic content of 
the blade response; deflections, slopes, moments, and shears; and 
of the control system response. The consideration of a deformed 
free-wake requires an iterative procedure to be applied to the 
fundamental analysis since the rotor system will respond to a 
downwash distribution resulting in a modified downwash distribu- 
tion on the rotor system. Wake-induced velocity influence co- 
efficients and the initial nondimensionalized bound circulation 
values necessary to determine the initial downwash distribution are 
obtained from use of an existing free-wake analysis. 

A digital computer program based upon the analysis developed 
was executed for a swashplate controlled, articulated rotor 
system in steady-state forward and maneuver flight configurations 
including a deformed free-wake to ascertain the viability of this 
theoretical analysis. The dynamic response results obtained from 
these calculations indicate that the coupled rotor system dynamic 
response analysis provides an accurate and improved simulation of 
a helicopter rotor system. 


*PART II - Computer Program Listing is NASA CR-2453. 



INTRODUCTION 


In recent years the design of helicopter and V/STOL rotor 
systems has required a more realistic representation of the 
dynamic and aeroelastic behavior of rotor systems. The pre- 
diction of rotor-induced aerodynamic flows has recently under- 
gone extensive refinement, primarily due to the need for im- 
proved predictions and to the development of large-capacity 
digital computer systems. The desire for increased performance 
and other considerations have stimulated the development of im- 
proved analytical tools for use in the design of helicopters. In 
spite of the technological advances which have been accomplished, 
helicopters continue to be limited in their operational envelope 
often as a result of severe oscillatory airloads and rotor system 
response. There still exists a need for a better understanding 
of the causative physical phenomena for the behavior of fully- 
coupled rotor systems. 

The analysis which was developed and the computer program 
which resulted are a part of an increased effort to develop 
better theoretical methods for the study of helicopter rotor 
systems. Improved high-harmonic airloads determination alone 
is of limited usefulness in blade response analysis unless 
compatible blade models are employed in the dynamic response 
computations. Inter-blade coupling due to the control system 
and aerodynamic coupling can have significant effects on the 
torsional dynamic response behavior of a rotating helicopter 
blade. In order that the behavior of helicopter rotor systems 
can be more accurately predicted, a theoretical analysis was 
developed which includes the effects of; high harmonic airloads, 
anisotropic control systems, interharmonic blade coupling, and 
coupling of blade motions with a deformed f ree-wake. 


SYMBOLS* 


a speed of sound 

th 

a m length of rigid rocker arm attached to m— — 

blade, positive if rocker arm aft of 
quarter chord, m 


*A11 units given in SI units but equivalent English units may 
be used if applied to all definitions of program variables. 
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[A], 


(i) 

m 





[B] (l) , [B] (l) 
m ' 1 1 m 


[B. ] (l) 

k,n J m 



,n^m 


{E 


(i) 


k,n m 



c 


j 


real time and Laplace transformed blade 
transfer matrix associated with mass and 

inertia effects for the i— station on the 

m— blade 

real part of a Np harmonic state variable 

th 

inboard of the i — station 

til 

rigid offset of j — linear spring swash- 
plate attachment from neutral axis of ring 
(positive if offset is toward the center 
of ring) 

real time and Laplace transformed blade 
transfer matrix associated with a bend or 
finite twist of a blade section 

transfer matrix relating contribution to 
the k-shifted state vector harmonic coef- 
ficients at the inboard end of section i 
of blade m from the n-shifted blade tip un- 
knowns harmonic coefficients, where, for 
example, a k shift denotes a harmonic 
-k/rev. relative to the Np/rev. harmonic 
(this double subscript notation denotes 
the same form of relationship in other 
arrays) 

portion of final set of rotor and swash- 

th 

plate equations associated with the m-— - 
blade 

matrix relating contribution to state 

th th 

vector for the m— - blade at i— — blade 
station from the control torque applied to 
blade through rocker arm 

imaginary part of a Np harmonic state vari- 

i, 

able inboard of the i— — station 

X. U 

linear damper strength of j— - cyclic 
spring-damper unit supporting the swash- 
plate, N-sec/m 
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'♦j 


Ci) 

m 


(c, '} 


(i) 


k,n m 


(£, ] 
k m 


[CJ 


(i) 

m 




(i) 

m 



n 


(i) 

m 



til 

lateral angular damper strength of j — 
cyclic spring-damper unit supporting the 
swashplate, N-sec/m 

longitudinal angular damper strength of 

+* Vi 

j — cyclic spring-damper unit supporting 
the swashplate, N-sec/m 


tii 

defines the blade chord at the i— blade 
station of the m-— blade 

matrix relating contribution to state 

vector at i— — blade station from the dis- 
continuity in torsional deflection at the 
pitch bearings 


portion of final set of rotor and swash- 
plate equation coefficients associated 

til 

with the coupling of the m-— blade with 
the swashplate response, defined by 
Equation (150) 

blade transfer matrix associated with a 
lumped aerodynamic station, accounts for 
periodically varying aerodynamic terms 
involving d/dt 

harmonic blade transfer matrix associ- 
ated with lumped aerodynamic station 
obtained by Fourier analysis of previous 
variable 


total aerodynamic blade transfer matrix for 

+- Vi 

a transfer across the i—~ station on the 
m~ blade 

lift coefficient 

lift coefficient limit for special stall 
model 


Cp lift curve slope at 0 radians angle of 

a' attack 


c^ drag coefficient 
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two-dimensional stall limit moment 


m 2D, 

stalled 


C 


Gx 


C 


Gy 


c 


c 


C' 


d 


m 


^ ( x j / i-) » d ( 6 , 1 ) 



} 


(i) 

m 


(d) 


(i) 

m 


(i) 
n' m' 




(i) 
n' m 


[U_ v J m 

n-K,n m 


damper strength of collective spring-damper 
unit supporting gyroscope ,N-sec/m 

lateral damping coefficient of gyroscope 
foundation support, N-m-sec/rad 

longitudinal damping coefficient of gyro- 
scope foundation support, N-m-sec/rad 

damper strength of collective (foundation) 
spring-damper unit supporting swashplate, 
N-sec/m 

average of longitudinal and lateral gyro- 
scope damping coefficient, N-m-sec/rad 

lateral gyroscopic damping minus the longi- 
tudinal damping coefficient /2, N-m-sec/rad 

fch 

rigid offset of m— - control rod attachment 
point from neutral axis of ring (positive 
if offset is outward from center of ring) , 
m 

definitions involving swashplate deflections 
to reduce length of equations 

matrix relating contribution to state vector 

th ... 

at i— — blade station from the discontinuity 

in flap angle at the flap hinge 

forcing function due to aerodynamics at an 
aerodynamic station 

real time and Laplace transformed harmonic 
aerodynamic forcing function obtained by 
Fourier analysis of the previous variable 

portion of final set of rotor and swashplate/ 
gyroscope equation coefficients associated 
with coupling of the swashplate or gyroscope 

response to that of the m— - blade, defined 
by Equation (148) 


5 



ID] 


(i) 

m 


] (l) [D ] 
n m, 1 n J 


(i) 

m 


[E] ( j } [E] 


e 

(i) 

m 


(EI y ) (i) 
(El ) (l) 

z 


(EJ Z ) 


blade transfer matrix of terms not multiplied 
by a time derivative for an aerodynamic 
station 

real time and Laplace transformed harmonic 
transfer matrix based on previous variable 

blade sweep hinge radial offset, m 

real time and Laplace transformed blade 
transfer matrix associated with a massless 
elastic section 

til 

local edqewise bending rigidity of i— 

2 

blade sect ion, N-m 

th 

local f lapwise bending rigidity of the i— 

2 

blade section, N-m 

2 

bending stiffness of swashplate, N-m 


F (a) 


{F°} ^ 
m 

{F 1 } ^ 
m 


variable involving swashplate stiffness to 
reduce equation length as defined after 
Equation (37) 

real time and Laplace transformed forcing 
function due to the mass and inertial effects, 
k subscript is added for shifted harmonics 


steady forcing function due to mass and 
inertia effects 


amplitude of +l/rev forcing function due to 
mass and inertia effects 


(F _1 } ^ 


amplitude of -1/rev forcing function due to 
mass and inertia effects 


{F^.}* total forcing function effect in final 

solution as defined by Equation (153) 

total swashplate forcing function, always 
zero 


{Fr^}* 



{F Pv } r 


totax blade forcing function in final solu- 
tion representation and as defined in 
Equation (154) 


gravity constant, m/sec' 


local torsional rigidity of i— blade 
section, N-m 2 


<GJ t ) 


tH. } 
k m 


torsional stiffness of swashplate, N-m 

real time and Laplace transform version 
of total section transfer matrix associated 

with the i— section of the m^ blade 

radial distance ( positive outward) of a 
blade mass point from the origin of the 
reference coordinate system at the rotor 
hub when 4> and 0 are zero in value , m 

distance of neutral axis ahead of pitch 
axis assumed to be at quarter chord, m 

total forcing function at the i~— station 

A. 

of the m— blade due to all outboard external 
forces and moments 


imaginary number, v r -T~ 


L dG'^dG' k dG 


r t r 
x f °f ,K f 


orthogonal unit vectors for gyroscope in a 
perturbed system 

orthogonal unit vectors for the fixed; shaft, 
swashplate, and gyroscope systems 


inertia dyadic of the gyroscope 


*t*" 

i 1 3 t k 
r J r r 


orthogonal unit vectors for the rotating? 
shaft, swashplate, and gyroscope systems 


y(i) q-(i) 
x rb '-'rb /K rb 


orthogonal unit vectors in local blade 

. th 

coordinate system at the i — station 
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I Q lateral moment of inertia of gyroscope, 

kg-m 2 

I„ longitudinal moment of inertia of gyroscope 

G Y ? 

kg-in 

I Qz polar moment of inertia of gyroscope, 

kg-m 2 



K 

K 


K' 



K b,k 


KC a 

^6 

a 



polar mass moment of inertia of a section 

2 

about the neutral axis, kg-m 

flapwise mass moment of inertia of a 
2 

section, kg-m 

edgewise mass moment of inertia of a 
section, kg-m 2 

stiffness of collective spring-damper unit 
supporting swashplate or gyroscope, N/m 

average cyclic stiffness of the gyroscope, 
N-m/rad 

lateral minus longitudinal cyclic stiffness 
divided by 2, N-m/rad 

variable to reduce equation length as 
defined after Equations (17.) and (18) 

variable to reduce equation length as 
defined after Equation (157) 

variable to reduce equation length as 
defined after Equation (67) 

variable to reduce equation length as 
defined after Equation (69) 

variable to reduce equation length as 
defined after Equation (69) 



k GX lateral angular spring stiffness of gyro- 

scope support, N-m/rad 

k QY longitudinal angular spring stiffness of 

gyroscope support, N-m/rad 

.fch 

k. stiffness of j — cyclic spring-damper unit 

3 supporting swashplate or gyroscope, N/m 

ttl 

k Q lateral angular spring stiffness of j — 

j swashplate or gyroscope support, N-m/rad 

kg longitudinal angular spring stiffness of 

j— swashplate or gyroscope support, N-m/rad 

k stiffness of the control rod attaching the 

m th 

m— - blade to the swashplate or gyroscope , N/m 

k stiffness of torsion spring, N-m/rad 

k stiffness of flapwise moment spring, N-m/rad 

k stiffness of chordwise moment spring, N-m/rad 

z 

( i ) tli 

L ' } length of the aerodynamic field for the i— - 

section if it has aerodynamics ,m 

£ ^ length of massless elastic section, m 

L ^ lift acting on the i— section for azimuth- 

al location denoted by k if aerodynamic 
section, N 

operator indicating Laplace transform with 
respect to time t 

m mass of swashplate, kg 

m^ blade section mass, kg 

moment acting on gyroscope from blade, N-m 
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M 

G 

M y ,M z ,T 

<V 


Mass of gyroscope, kg 

f lapwise bending moment, edgewise bending 
moment and pitching moment in a rotor 
blade, N-m 

azimuthal harmonic of bending moment in 
swashplate, N-m 


M 


sp 

"z 

N. 


N c . t . ,N f ea' 


N 


flap 


N 


e . s . 


N , 


N 


ff 


bending moment in swashplate, N-m 
number of rotor blades 

blade station immediately inboard of which 
the rocker arm, pitch bearings, and flap 
hinge occurs, respectively 

number of concentrated supports below swash- 
plate or gyroscope 

maximum number of harmonics included for 
interharmonic coupling 

maximum number of blade loading harmonics to 
be used 


N 

max 

N p ,NP 

NS 


P 


(i) 




{P„}„ 

n ro 


highest azimuthal harmonic of swashplate 
deformation accounted for 

harmonic response of interest 


number of stations or segments into which a 
blade is divided 


centrifugal force components acting in the 
local spanwise direction, N 

real time and Laplace transformed counter- 

1.1. 

part of force in control rod of the m~— 
blade, N 

column of blade tip unknowns and blade 
discontinuities 

variable involving swashplate stiffnesses 
defined to reduce size of equations 



Q(e,t) 

Q £ ( t ) 


R 


m 


m <*> m <*> 



} 



distribution o± forces acting on swashplate 
in a downward sense, N/m 

- fourier harmonic of the loading 
Q (6 ,t) , N/m 

column containing dependent variables of all 
blades plus swashplate or gyroscope 

blade radius from axis of rotation, m 

position vector of blade attachment point 
to the gyroscope 

real time and Laplace transformed blade 
transfer matrix associated with a rigid 
massless segment-may be used for an offset 
in neutral axis 

column containing harmonics of swashplate 
deflection or rigid body response of 
gyroscope 


amplitude of a Np/rev variable of interest 
at the i-— - section 


R 


R 


sp 


3 




s 


S 


a 



] 


radius of gyroscope, rn 
radius of swashplate, m 

row vectors having 12 elements, all of which 
are zero except for a value of unity in the 
third, fourth, and eleventh positions, 
respectively 

complex number which is the Laplace trans- 
form variable corresponding to time t, 
rad/sec 

variable to reduce equation size defined 
following Equation (37) 

portion of the final set of rotor and swash- 
plate equation coefficients associated with 
the swashplate impedance, defined by 
Equation (147) 
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(S } ^ , { S } ^ real time and Laplace transformed version 

n 1 " 11 of the state vector at station i of the 

til 

m— blade, subscript n meaning the Laplace 
transform variable s is replaced by s-infi ,- 
added to superscript denotes vector at 
inboard end of section where + denotes 
outboard 

{S Laplace transformed state vector for n 

n m th 

shift and m— - blade consisting of all 

twelve state variables 

Laplace transformed modified state vector 

4* V* 

for n shift and m— - blade consisting of 
six non-zero tip state variables 

real time and Laplace transformed version 
of blade transfer matrix associated with 
concentrated spring-damper unit (lead-lag 
spring-damper unit) 

t time , sec 

[T^ ] contains coefficients for final set of 

' n equations governing swashplate or gyroscope 

response, blade tip unknowns and blade 
discontinuities, defined by Equation (146) 

azimuthal harmonic of swashplate torque, N-m 

T S P torque in the swashplate, N-m 

u 0 (t)fU Q (s) real time and Laplace transformed displace- 
ment downward of swashplate base, m 

u ,u ,u perturbation deflections of blade neutral 

axis in radial, edgewise, and flapwise 
senses , respectively , m 

v ^ induced velocity at i— section if uniform 

azimuthal ly, m/sec 


[SK] , ISK] { t ) 
m m 
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v(e,t) 



v 


k 


(i) 


v,(t) 


VV N 


displacement downward of swashplate as a 
function of azimuth and time (referred to 
fixed frame) , m 

induced downwash (inflow) velocity at the 

i— blade section and j— azimuthal station, 
m/sec 

tin 

freestream velocity at the i— section and 
th 

k— - azimuthal location, m/sec 
azimuthal harmonic of v(6, t) 

edgewise shear force, flapwise shear force, 
and radial force in blade, respectively, N 


<V» 


azimuthal harmonic of V 


sp 

y ' 


N 


V S P 

Y 


transverse shear force in swashplate, N 


W' 11 ' 


V+k> 


(i) 


V ,N 

Z ' p 


chordwise and flapwise shear forces acting 
at inboard end of section (1) at the 
azimuthal position t^, N 

Np/rev harmonic of the flapwise shear force 
at the inboard end of section (i) , N 


V 


z 


(i) 


(t) 


real-time flapwise shear force at inboard 
end of section (i) , N 



gyroscope impedances associated with 
asymmetry of support configuration 


V Q forward speed of helicopter, m/sec 


w(<(),t) displacement of swashplate as a function of 

azimuth and time (referred to rotating 
frame) , m 

w £ (t),w^(s) azimuthal harmonic of w(<j>,t), real-time and 

Laplace transformed versions for swash- 
plate, respectively 
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azimuthal harmonic of w(<j>,t,), real time 
and Laplace transformed version for gyro- 
scope, respectively 

rectangular Cartesian coordinate nonrotating 
shaft system 

perturbed gyroscope coordinate system 

rectangular Cartesian coordinate system 
rotating about z^ axis at n radians/second 

local right-handed coordinate system attach- 
ed to blade, Y ^ directed along chord, 3 ^ 

normal to blade 

swashplate impedances associated with 
asymmetry of support configuration 

gyroscope impedances associated with 
asymmetry of support configuration 





[a] 


“is 

a lc 


complex function characterizing rocker arm 
length and control rod retardation time 

swashplate impedances associated with 
symmetric properties of swashplate and its 
support configuration 

\ 

gyroscope impedances associated with 
symmetry of support configuration 

total blade section angle of attack at the 
i— section and k— azimuthal position, rad 

boundary condition array on blade root dis- 
placements and slopes 

amplitude of cyclic input, rad 

lateral cyclic angle, rad 

longitudinal cyclic angle, rad 
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a s 

shaft tilt angle, rad (ir/2 for hover) 


3 

cyclic phase angle, tan X ( — — ) f rad 

" a ls 


3 ( 8 r t) 

applied torque per unit length acting on 
ring, N-m/m 


*»<*> 

Fourier harmonic of g(e,t), N-m/m 


6 U) 

B Np 

phase angle of state variable as shown in 
Equation (167) , rad 


r 

nondimensionalized circulation 


r (i) 

! k 

i, 4- 

circulation at the i— section and k— 
azimuthal position, m/sec 



kronecker delta, equals unity if i= j , 
equals zero if i/j 


6( ) 

Delta function 

{"a } c * ^ * 
J m 

{ij* iap 

m 

- ■ 

discontinuity columns defined by Equations 
(91), (87) and (89), respectively 

— *c.t. - ,ffea 

^ A n ^m ' ^ A n^m 

a i* flap 

1 n vm 

discontinuity single element arrays defined 
by Equation (95) with Laplace transform 
variable shifted by -ing 

(at) . 
nr 

(Ad) ) 
y y m 

(Ait ) , 

y x'm' 

discontinuities in blade torque, torsional 

til 

deflection and flap angle of m — blade 


e ry 

distance of mass center of gravity ahead of 
elastic axis, m 


6 

azimuthal independent coordinate for swash- 
plate unknowns referred to fixed frame of 
reference, rad 


0(0, t) 

local bending slope of elastic ring repre- 
sentation of swashplate, rad 
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collective pitch, rad 


6 c 

0^(t) Fourier harmonic of 0 (9 , t) , rad 

{ Xj, n > m initialization array with xxx superscript 

' equal; c.t. , fea , and flap ; correspond- 

ing to control torque, feathering angle, 
and flap angle, respectively, as defined 
by Equation (96) for harmonic k and shifted 
by -inft 

12 by 6 matrix utilized to eliminate zero 
value tip variables 

\i mass per unit arc length of swashplate, kg/m 

p air mass density, kg/m 3 

it Pi, 3.1415927 

[a] wake-induced velocity influence coefficient 

matrix 

a(e,t) applied moment per unit length acting on 

swashplate, N-m/m 

a (t) Fourier harmonic of o(0,t), N-m/m 

X/ 

damper time constant associated with spring- 
damper unit representation, of control rod 
t h 

connecting m— - blade to swashplate, sec 

$ , 0 , ¥ finite angles defining orientation of blade 

surface at an arbitrary radial point, $ 
corresponds to forward sweep, 0 corresponds 
to downward coning, V corresponds to nose- 
up pre-twist and/or collective angle of 
attack, rad 

azimuthal independent coordinate for swash- 
plate motion referred to rotating frame, of 
reference, rad 

$(e,t) swashplate local twist, rad 

(t) Fourier harmonic of 4>(6,t), rad 

X, 
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<f> m phase angle of m- — blade ( 2 tt (m-1) /N^ for 

equally spaced blades) , rad 

<f> m p perturbation on ^ or blade wrap-up angle 

on the m— blade, rad 


^x'^v'^z perturbation twist, flapwise bending slope, 

y and edgewise bending slope respectively, 

rad 



gyroscope perturbation variables, rad 

til 

azimuthal angle of j — spring-damper unit 
supporting swashplate , rad 


^ ^ blade azimuth angle, measured from down- 
stream position plus tt/2 in direction of 
rotation, rad 


ft operating speed of rotor system, rad/sec 

io o vortex interaction frequency 

Subscripts and Superscripts: 

tt 

i or (i) superscript denoting i— blade station 

counting inboard from tip to root, has 
successive values from 1 to NS 


j index 


£,n,p,q,k 

m 


+ - 


c. t. ,fea ,f lap 


( ) 


indices used for harmonic number 

th 

index indicating m — blade 

employed as superscripts to indicate the 
outboard and inboard ends of a blade seqment 
or radial station, respectively 

employed as superscripts to indicate a 
matrix product or state vector evaluated at 
the control torque, pitch bearings, and flap 
hinges, respectively 

as in (r 3 ) to indicate a row vector 


[ ] indicates a matrix 

{ } indicates a column vector 
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HELICOPTER ROTOR SYSTEM MODEL 


The analysis to be presented is dependent upon the defi- 
nition of the coupled helicopter rotor system model. Basically, 
the rotor system configuration is comprised of a rotor consist- 
ing of an arbitrary number of flexible blades, arbitrarily 
located azimuthally , with either a flexible swashplate or gyro 
control system. To provide model flexibility such that the 
analysis will be applicable to various rotor systems, the re- 
presentation of each of these major components consists of several 
relevant characteristics. 

The model of the blades comprising the rotor allows for the 
inclusion of all blade characteristics believed to be significant. 
For example, the blade model allows: 

(a) arbitrary location of the blade elastic axis and 
aerodynamic center , 

(b) arbitrary location of center of mass relative to blade 
elastic axis, 

(c) arbitrary blade precone and presweep, 

(d) arbitrary variable twist distribution, 

(e) arbitrary mass and inertia distributions, 

(f) arbitrary stiffness distributions, 

(g) gyroscopic forces 

(h) aerodynamic loading including damping, 

(i) flapwise, chordwise, and torsional blade root restraints 
ranging from rigid to fully articulated and, also, 
blade root damping. 

A lumped parameter approach to the modeling of the blades has been 
adopted as a suitable means of representation of the blade distri- 
buted properties. 

The main component of the swashplate control system model is 
represented by a flexible ring consisting of two parts which is 
allowed to translate only along the shaft axis but may rotate 
about two mutually orthogonal axes perpendicular to the shaft 
axis. The upper portion of the ring is also allowed to rotate 
about the shaft axis since it is attached to blade pitch arms by 
control rods having both stiffness and damping characteristics. 

The lower portion of this ring is supported by an arbitrary number 
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of supports which have stiffness and damping characteristics and 
are attached to a ring allowed only translational motion in the 
shaft direction (collective motion) . The support for this final 
ring has an arbitrary effective stiffness and damping value. 

This final support is necessary since the collective stiffness 
can be much lower than the cyclic stiffness of the control 
system. Thus, by variation of stiffness and damping characteri- 
stics of the supports involved any degree of anisotropic support 
of the swashplate control system can be obtained. Additional 
capabilities included in the swashplate control system model are 
presented in a subsequent section. 

The gyro control system model considers a rigid gyroscope 
that is rotating about its center of gravity which is allowed to 
translate in the shaft axis direction. The gyroscope is supported 
in such a manner that arbitrary lateral, longitudinal, and col- 
lective stiffness and damping values can be associated with the 
support system. The gyroscope is attached to the blade pitch arms 
by control rods having stiffness and damping characteristics. 
Additional information on this control system configuration is 
presented in a subsequent section. 

The general rotor blade coupling to a control system is 
shown schematically in Figure 1, for a swashplate control system. 
Also depicted in this figure are aerodynamic load points and 
their attached trailing vortices representation and the relation- 
ship between the rotating shaft coordinate system (r-subscripted) 
and the nonrotating or fixed shaft coordinate system (f -sub- 
scripted) which will be involved in the analysis to be presented. 


GENERAL SOLUTION METHOD 


In contrast to computation of blade response by use of normal 
modes of free vibration as generalized coordinates , the method of 
analysis that is used is based, upon a direct inverse of a dynamic 
matrix. This method eliminates the need for normal mode fre- 
quency and mode shape information in determination of the dynamic 
response of a helicopter rotor system and also eliminates various 
engineering judgments such as the number and types of normal modes 
which must be included to yield proper results . The development 
of the dynamic matrix and its associated forcing function and 
subsequent solution by taking a direct inverse is straightforward. 
The basic method of solution can be easily adapted to provide 
coupling of a deformed free-wake to the blade motions by an iter- 
ative procedure. 

The analysis for representing blade behavior in the dynamic 
matrix generation is accomplished by application of a transfer 
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Figure 1. Coupled blade-swashplate representation 



matrix approach in the lumped parameter blade model. Basically, 
the usual transfer matrix procedure consists of the following: 

(a) a modeling procedure in which the distributed blade 
properties are converted to a lumped parameter form 
representing the blade as consisting of uniform 
massless elastic beam sections, in tension due to 
centrifugal loading with point mass , inertias , and 
aerodynamics located at the ends of the massless 
lengths , 

(b) a process in which section transfer matrices associated 
with successive sections of a blade and relating shears, 
moments, slopes, and deflections at each end of the 
section are multiplied together to form the transfer or 
associated matrix for the total blade. 

The transfer matrix procedure used in this analysis is a 
modification of that outlined above. The usual transfer matrix 
technique involves arrays of real numbers since the governing 
equations from which the transfer matrices are developed are 
ordinary linear differential equations with constant coefficients 
and have no damping or gyroscopic terms. For natural frequency 
prediction, the second time derivative may be replaced by the 
negative of the natural frequency squared. However, for an 
n-bladed rotor system which is subject to aerodynamic, gyroscopic, 
and damping forces, the governing differential equations have 
periodically varying coefficients such that it is no longer 
possible to simply replace the differential operator by a natural 
frequency term. By applying the Laplace transform to all vari- 
ables, including periodically varying coefficients and different- 
ial operators, in the real-time transfer matrices, the transfer 
matrices are altered to arrays of complex numbers. The process 
of combining the complex transfer matrices to form the associated 
matrix module representing the rotor blades is straightforward 
but differs considerably from the common transfer matrix procedure 
because of harmonic coupling introduced by the aerodynamic time- 
varying coefficients. 

Application of the modified transfer matrix procedure across 
the blade sections from blade tip to blade root yields the de- 
finition of the shears, moments, slopes, and deflections at the 
blade root in terms of the corresponding unknown tip values. 
Boundary conditions that are used in the generation- of the final 
blade governing equations require that the shears and moments at 
the blade tip be zero and that the slopes and deflections at the 
blade root are zero. 

The anisotropically supported swashplate or gyroscope control 
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system is represented by analytical expressions for the control 
system motion obtained from application of the differential 
equations of continuity and equilibrium to the control system 
model. Since the blade and control system reference coordinate 
systems differ, the resulting equations of motion which include 
coupling to the blades thru the control rod forces are trans- 
formed to the blade reference coordinate system to provide com- 
patibility of force application. Also, the control system equa- 
tions of motion are Laplace transformed to replace differential 
operators with a complex number representation compatible with 
the blade analysis. 

The control system equations and the blade equations re- 
sulting from application of blade boundary conditions are used 
to construct a dynamic matrix and a forcing function column. 

The solution for the unknown tip quantities and the control system 
unknowns is obtained by multiplication of the forcing column by 
the inverse of the dynamic matrix. From a knowledge of the blade 
tip unknowns the transfer matrix procedure is repeated across the 
blade sections to obtain the radial distributions in complex 
number form representing both amplitude and azimuthal phasing 
for shears, moments, slopes, and deflections. This method of 
solution obtains the blade response by harmonics of the rotation- 
al speed where the value for the Laplace transform variable de- 
termines for which harmonic the analysis is being solved. The 
total real-time radial and azimuthal distributions of a dynamic 
response variable is obtained by a proper summation of harmonic 
components of that variable. 

A general representation of the solution method without free- 
wake coupling included is presented in flow chart form in 
Figure 2. The notation," intermediate terms", denotes terms that 
are not a function of the Laplace transform variable and there- 
fore, need not be recalculated for different harmonics of dynamic 
response. Downwash velocities may be taken as uniform, radially 
varying, or both radial and azimuthally varying in this solution 
form in which the blade response includes downwash velocity 
effects on the aerodynamic terms but coupling of the free-wake 
and blade motion is not allowed. 

Deformed free-wake effects on the helicopter rotor system 
dynamic response are included in the analysis by modifying the 
solution method to include an iterative procedure which couples 
the free-wake to the blade motions. The iterative procedure 
developed requires the execution of a free-wake analysis program 
(wake geometry program) developed for reference 1 but a similar 
free-wake analysis program could be substituted if it provides 
the necessary information. The wake geometry program with 
flapping and cyclic blade positions defined from either an airloads 
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calculation or experimental values uses a real-time analysis to 
generate a fully deformed free-wake corresponding to the blade 
positions. In addition, the wake-induced velocity influence co- 
efficient matrix, [a] , and the nondimensionalized bound circula- 
tion distribution, { r } , are determined and available for use in 
the dynamic response iteration procedure. 

The a ' s and r ' s from the wake geometry program are used to 
determine the initial wake-induced velocity (downwash) at the 
blade radial and azimuthal loading points (i. e. [a] { T } « (v) ) . 

This initial downwash distribution is introduced into the aero- 
dynamics analysis in which associated angles of attack and 
resultant velocities are used to determine the harmonics of blade 
aerodynamic loadings. With these loadings, the harmonics of 
dynamic response including control system effects from steady 
response up to the highest harmonic of response to be included in 
the wake iteration procedure are determined. But the dynamic 
response of the system includes blade motions which alter blade 
aerodynamic loading. From consideration of the total shear re- 
sponses acting on both sides of the blade aerodynamic load points, 
radially and azimuthally located, a lift distribution can be 
obtained which provides a new bound circulation distribution, 

{ r } . A new downwash distribution is obtained by use of the 
initial wake-induced velocity coefficient matrix and the new 
{ r } . This new set of downwash values which includes coupling of 
the free-wake to the blade motions is then inserted into the 
aerodynamics analysis and the procedure repeated. 

The procedure involved from one set of downwash values to 
the obtaining of a modified set of downwash values (an iteration) 
can be repeated until a specified maximum number of iterations 
between blade circulations and blade motions are performed or 
the values of blade circulations for successive iterations agree 
to within a certain specified value (convergence criteria) . 

When either of these conditions occur, the harmonic content, in 
complex form, of the blade moments, shears, slopes, and deflec- 
tions along the blade and of the swashplate motion can be obtained. 

A general representation of the solution method with free- 
wake coupling to blade motions included is depicted in flow chart 
form in Figure 3 where the segment with the notation, "dynamic 
response segment"/ represents operations which are part of the 
basic solution method and is labeled on Figure 2. Aerodynamic 
analysis must be included in the free-wake iteration solution 
method. 


24 




25 


















ANALYSIS 


Derivation of the Governing Equations for an 
Aniso tropically Supported Gyroscope 

In order to allow flexibility in the choice of control systems 
which may be utilized in the analysis representing a coupled heli- 
copter rotor system, the analysis for a form of a gyroscopic 
control system is developed. This analysis provides the necessary 
representation for the gyroscope control system behavior and for 
the coupling which occurs between this control system and the 
blade behavior. 

Gyroscope control system model .- The model chosen for the 
development of analysis representing a gyroscope control system 
consists of a rigid gyroscope, which may be nonsymmetrical , that 
is supported by a set of spring-damper units to allow anisotropic 
support. In addition, the gyroscope is attached to the rotor 
blades by flexible control rods having both stiffness and damping 
characteristics. The complexity of the analysis is restricted by 
assuming the gyroscope to be allowed only translational motion 
of its center of gravity in the shaft direction and rotational 
motion about three mutually orthogonal axes located at the center 
of gravity. Analysis based on these assumptions should represent 
the major portion of the gyroscope control system effects on the 
blade dynamic response behavior. 

Coordinate systems invo lved.- A development of analysis to 
represent the' gyroscopic control "system behavior requires the de- 

finition of three mutually orthogonal coordinate systems. These 
three coordinate systems consist of a fixed (stationary) system, 
a rotating system, and a perturbation (local) gyroscope coordinate 
system. 


The fixed gyroscope coordinate system is located in space 
such that the gyroscope center of gravity location in an unper- 
turbed situation (point 0) coincides with its origin as shown in 

Figure 4 where the fixed system unit vectors T^, and corres- 

pond to the x f , v f , and z f axes, respectively. This fixed coordi- 
nate system has its axes oriented in the same directions as the 
shaft fixed coordinate system which will be used in later sections. 
These two fixed systems are different in that their origin loca- 
tions may occur at different positions in the shaft direction. 

The rotating gyroscope coordinate system has the same origin 
and z axis as the fixed coordinate system but its x and y axes 
are rotating at a rotational speed of ft relative to the fixed 
system x and y axes. The rotating coordinate system axes and unit 
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vectors are denoted by r subscripts as shown m Figure 4. Also 
shown in Figure 4 is a top view of the coincident x-y planes of 
the fixed and rotating gyroscope coordinate systems from which 
the rotating coordinate system can be seen to be related to the 
fixed coordinate system by the coordinate transformation 
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where t represents time and fit represents the instantaneous angle 
between the two systems. The rotating gyroscope coordinate 
system also exhibits the same relationship to the shaft rotating 
coordinate system which will be used in later sections as was 
noted for the fixed gyroscope system versus the fixed shaft 
system. 


The local gyroscope coordinate system, dG subscripts, has its 
origin at the center of gravity of the gyroscope (point G) and 
translates as the gyroscope translates. The z^ G axis always co- 


incides with the principal axis of inertia of the gyroscope that 
is perpendicular to the plane of the gyroscope. The other two 
mutually orthogonal axes lie in the plane of the gyroscope such 
that if the plane of the gyroscope is parallel to the rotating 
coordinate system the x and y axes of the local coordinate system 
would be in the same direction as the x and y axes of the rotating 
coordinate system. This coordinate system is presented in general 
concept with the other two coordinate systems in Figure 4. 

The orientation of the local gyroscope coordinate system 


z^ G ) relative to the rotating coordinate system (x^, y r , z^) is 
obtained by performing three successive orthogonal rotations; $ 

<f> , and 4> ; on the rotating coordinate system as depicted in 

V X 

2 x r 


z 


r 


f 


Figure 5. In terms of the unit vectors of the two coordinate 
systems and the perturbation rotations relating these coordinate 
systems, the relation of the local gyroscope system to the rotating 
gyroscope system assuming small angle approximations is expressed 
by the coordinate transformation 
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Neglecting second order terms 


the inverse transformation is 
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Real-time analytical derivation. - The real-time analysis for 
representation of the gyroscope control system requires the de- 
velopment of the equations of motion of the gyroscope. These 
equations of motion are developed considering the perturbation 
rotation related equations to be independent of the translational 
equation of motion. In regard to the analysis presented in this 
section; notation of the form, a, is used to represent the first 
time derivative of the variable, a, and notation of the form, 

a, is used to represent the second time derivative of the variable, 
a. In addition, the time derivative of a vector quantity includes 
a contribution due to the vector product of the rotational vector 
and the vector of interest. 

The rotational velocity (angular velocity) vector of the ro- 
tating coordinate system can be represented by 

n = n . 

The angular velocity of the gyroscope is given by 


ft = ftk,. + $ l+l j + k 

gyro f y x r J r z r 


( 5 ) 


referred to the x , y r , rotating coordinate system or 


fi gyro " (4> x " *y fi) i dG + (<l> y + ^x n) j dG + {<f> z + ^ )k dG (6) 

VJL r ; r J r r r 
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as referred to the body local system of axes with the components 
being along the instantaneous direction of the axes. 


The inertia dyadic of the gyroscope at its center of gravity 
G, referred to the axes x^ G , y^ G , z^ g is 


I Gx X 4G ^dG + I Gy ^dG ?dG + J Gz k dG k dG 


(7) 


The angular momentum vector of the gyroscope at its center 
of gravity (principal axes) i^ from Equations (6) and (7) 

L„ — I _ (4> ” fi) i + I_ (<|> + <j> ft)j,_ + I ( 0 , + <(> ) k,_ 

-G Gx x r Y y r dG Gy y r x r J dG Gz z r dG * 

( 8 ) 

Referring to the x^., y , z^, rotating coordinate system the angular 
momentum is given by 
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or rearranging 
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Then, it follows that 
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The moment acting on the gyroscope from the support system is 
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with K Gx , K q being the lateral and longitudinal stiffnesses of the 
gyroscope support system and C Gx , being the lateral and 

longitudinal damping coefficients of the gyroscope support system. 

The moment acting from the blade on the gyroscope is 

N, 


"b = ' l. r m x P m (t)k j 

m=l 


(13) 


or replacing r with its equivalent 
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(14) 


b 

M. = R y P Ct) I cos Cd> ~ 4> li — sin(d> ~ <j> )i ] 
d p m r m Y mp J r y m r mp r 

where is the number of blades , 

th. 

P m (t) is the force in the m— ■ control rod, 

Rp is the radius of the gyroscope, 

4> is the starting azimuthal position of each blade, and 
<j) m p is the perturbation on <j> m (blade wrap-up angle) . 

The governing equation of motion of the gyroscope for rota- 
tions is 

ic = M s + V <1S) 

The governing equation of motion of the gyroscope for vert- 
ical translation is 
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where M_ is the mass of the gyroscope, 

K is the collective stiffness, and 
C Q is the collective damping. 


The final real-time equation of motion of the gyroscope rota- 
tions is obtained in vector form by insertion of Equations (11) , 

(12) , and (14) into Equation (15) and replacing the sin(2fit) and 
cos(2st) by their exponential equivalents. Since Equation (15) 
must hold on a vector component basis, two expressions are obtained 
by considering only the terms multiplying the T and J , independ- 
ently. The resulting real-time equations of motion for gyroscope 
rotations are 
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+ [ aCI Gz - *Gx - V It + iK J ~ 1K -1 - a CJ *y r 
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where for convenience 

Kj = e 2aiQt (K' + b^C’i + C' |^)/2 . 

The final real-time governing equation of motion for gyroscope 
translation is obtained by rewriting Equation (16) in the form 
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Application of Laplace transforms .- The final real-time gyro- 
scope governing equations are Laplace transformed to eliminate the 
derivative type of representation and replace it with an alge- 
braic form in addition to providing compatibility with the blade 
equations which will be developed. In taking the Laplace trans- 
forms of the governing equations of motion quiescent boundary 
conditions are assumed and the following Laplace transformations 
required: 
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(d) L t [e lat FCt)_| = F(s-ia) 

(e) L t [e iat F(t)] = (s-ia) F(s-ia) 

where F(t) is a function of time. With these definitions, the 
Laplace transformed versions of the necessary equations can be 
obtained. 

The Laplace transformed versions of Equations (17) , (19) , 
and (18) are 
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where for convenience 

K b = (K' + bftC-'i + C's)/2 

The control rod forces P (s) couple the rotor blade motions 

to the gyroscope motions. The relations expressed by Equations 
(20) , (21) , and (22) are the objective of this section and will be 

used in a later section of this report. 


Derivation of the Governing Equations for an 
Anisotropically Supported Swashplate 

The form of control system utilized in the majority of heli- 
copter rotor systems is of a swashplate type. In this section the 
analysis for a general form of swashplate control system is de- 
veloped. This analysis provides the necessary representation for 
the swashplate control system behavior and for the coupling which 
occurs between this control system and the blade behavior. 

Swashplate control system model. - Even though the basic swash- 
plate model has been previously described there exist additional 
model capabilities. The general swashplate control system model as 
shown in Figure 6 consists of: 

(1) a flexible ring having uniform mass consisting of two 
parts which are allowed to translate only along the shaft 
axis but may rotate about the mutually orthogonal axes 
perpendicular to the shaft axis and where the upper por- 
tion is also rotating about the shaft axis, 

(2) a set of supports which have stiffness and damping 
characteristics which allow anisotropic supporting of the 
flexible ring and which may be attached to the ring an 
offset distance from the ring radius, 

(3) a single collective support which has stiffness and damp- 
ing characteristics in series with the previous group, 

(4) control rods with stiffness and damping characteristics 
connecting the blades to the flexible ring and which may 
be attached to the ring an offset distance from the 
ring radius , 

(5) a collective base plate which is allowed to only trans- 
late in the shaft axis direction, and 

(6) torsional spring-damper units which counteract the local 
bending slope and twist of the ring at each azimuthal 
location of the supports of (2) (shown on Figure 6 at 
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Figure 6 . Swashplate system model 


only one location) . 

Due to the generality of the configuration various swash- 
plate systems can be represented by- choice of the parameters 
used in the swashplate control system analysis. 

Coordinate systems involved. - Two mutually orthogonal co- 
ordinate systems are utilized in the development of the analysis 
representing the swashplate behavior. The first coordinate 
system is a fixed (non-rotating) coordinate system whose orienta- 
tion in space is identical to the shaft (hub) and gyroscope fixed 
coordinate system. However, the origins of these three fixed 
coordinate systems may be at different positions on the shaft axis 
(shaft centerline) . In these systems; the y-axis is perpendicular 
to the shaft axis and in a direction toward the nose of the heli- 
copter, the x-axis is perpendicular to the shaft axis and in a 
direction toward the helicopter's starboard side (advancing blade 
position) , and the z-axis is coincident with the shaft axis. The 
swashplate fixed coordinate system has been included in Figure 6. 

The second swashplate coordinate system is a rotating coordi- 
nate system whose origin and z-axis coincide with the swashplate 
fixed coordinate system but is rotating at a rotational speed, 

Q , such that the same relationship between the rotating and fixed 
swashplate coordinate systems exists as was depicted for the 
gyroscope coordinate systems in Figure 4. The swashplate, gyro- 
scope, and shaft rotating coordinate systems all have the same 
orientation as a function of time. 

Real-time analytical derivation. - The governing equations of 
motion of an elastic ring in terms of variables which are defined 
relative to the fixed swashplate coordinate system are used to 
obtain a governing differential equation for deflection behavior 
of the swashplate relative to the swashplate rotating coordinate 
system. The final form of the governing equation can not be ob- 
tained without the use of Laplace transformations. 

Fixed-frame governing equations involving loading functions : 
The governing differential equations for motion of an elastic 
ring involving variables which are a function of e, the azimuthal 
location of a point of interest on the ring relative to the 
x^-axis, and time, t, consist of the equilibrium, stress-strain, 

and strain-displacement equations: 
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V sp , M® p , and T sp are shear force, bending moment, and torque, 
y z 

respectively; 

v(0,t) is the local displacement of the swashplate; 

y is the mass per unit arc length of the swashplate ring; 

0(0,t) and $(0,t) are the local bending slope and twist 
angle of the ring; 

R g p is the swashplate radius; 

EJ is the bending stiffness of the ring; 

GJ t is the torsional stiffness of the ring; 

a(6,t) is the applied moment per unit length acting on 
the ring; 

0(0,t) is the applied torque per unit length acting on 
the ring; and 

Q(0,t) is the applied force per unit length acting on the 
ring . 


39 



The effects of concentrated moments or torques applied to the 
swashplate on the swashplate behavior are incorporated in terms 
appearing in Equations (24) and (25 } . The orientation of each 
of the variables used in Equation (23) through (28) are shown in 
Figure 7 . 



Figure 7 . Coordinates of swashplate variables 

The dependence of the governing swashplate partial differ- 
ential equations may be removed by going to a Fourier transform. 

Since each of the quantities v(0,t) , $(0,t), Q(d,t), T S & , and 

M S ^ are periodic with respect to 0 , each of these variables satis- 
z 

fies relationships of the form exemplified by v( 0 ,t) in the two 
Equations, (29) and (30) . 
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Thus, Equations (23) through (28) are multiplied by 
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the integration is performed over 0 from 0 to 2y . After using 
the conditions of periodicity and the definitions exemplified by 
Equation (30) above, the following set of ordinary differential 
equations with time as the single independent variable result. 


i £ -JLJL = 


R 


sp 


9 2 v. ft) 

* ~ 2 Q t (t> 

at 2 


(31) 


i£ 


<v t 


R 


sp 




(32) 


i£ 


R 


sp 


( Vi 


R 


sp 


- V (t) 


(33) 


i£ 


R sp 


*£ (t) + (M Z } £ 
■ + 


R 


sp 


EJ 


(34) 


v r (t) 

111 -r~ - 

sp 


(35) 


i£ 


V * 1 


R 


sp 


e £ (t > 

R sp 


GJ , 


(36) 


The same set of equations are also obtainable by direct substitu- 
tion of variable definitions based on the form of Equation (29) in- 
to Equations (23) through (28) and taking the equations correspond- 
ing to a specific value of £. 

Equations (31) through (36) can be combined to construct the 
following equation for v^ (t) in terms of the applied loading 
variable harmonics. 


3 v £ (t) EU) 

h + EML v ( t ) 

3t 2 R 3 £ 


Q,(t) 


i£ 


sp 


R 


a £ (t) + e £ (t) (37) 


sp 
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where for convenience 


and 


F(a) 



Ca 2 -1) 



2 




a 2 (GJ. + EJ ) 
t 2 

fa 2 GJ t + EJ z ) R sp 


When £ = 0 or ±1, F (£) and S are zero and swashplate stiffnesses 

A/ 

have no effect on v (t) since these are rigid body degrees of 

Jo 

freedom of the swashplate. During the manipulation of Equations 
(31) through (36) to obtain Equation (37) , additional equations 
can be obtained to represent the remaining swashplate slopes, 
moments, and shear variables in terms of the swashplate deflect- 
ion and applied loading functions. These equations are 


( v* = i Hr7v« 

sp 


(M ) „ = - zkl v, (t) + TP T- gg - $„ (t) 


z £ 


T, 


R P„ 
sp £ 

a ( £ 2 —i) 


GJ t P £ £ 


i£ R 


£ „ (t) 


2 V £ Ct) + EJ P P £ 

R sp P £ 


0,(t) 


R„ 


#„(t) = S„ v f , (t) + =ry — B»(t) 


£ £ 

i£ , 
r— 

sp 


EJ z GJ t p £ £ 


where for convenience 


(37a) 

(37b) 

(37c) 

(37d) 

(37e) 
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The variables defined by Equation (37) and Equations (37a) 
through (37e) and the loading functions involved in these express- 

th. 

ions correspond to the amplitude of the harmonic of the related 

time dependent variable at the azimuthal position on the swash- 
plate corresponding to the x f axis of the fixed coordinate system 

(6 = 0) . However, in order to couple the swashplate motion with 
the motion of the rotor blades , it is necessary to have both 
motions determined relative to either fixed coordinate systems or 
rotating coordinate systems. It has proven to be more convenient 
to use the swashplate deflections measured relative to a rotating 
reference frame corresponding to the blade rotating frame. 

Governing deflection equation in rotating frame involving 
fixed frame loading functions : By defining 

w(<£ ,t) = v(fit + <f>,t) (38) 

where fi is the operating speed of the rotor, the quantity w(4>,t) is 
equal to the downward vertical displacement of the swashplate at 
time, t, measured at an azimuth angle, <f> radians, in advance of the 
reference blade. This assumes that the reference blade makes an 
angle fit with the x-axis as shown in Figure 6. 

By exercising the Fourier representation as shown in Equation 

(29) 

oo 

w(<J>,t) = l w (t) e 1 ^ 

and 

oo , 

V(0t + ♦,« = l v,Ct> e liCot+ * ! '. 

Substitution of these expressions into Equation (38) and taking 
advantage of the characteristic that the resulting equation must 
hold for each harmonic 

, , , i£ (fit+<j>) .... i£<f) 

v £ (t) e y ' = w £ (t) e 

or in an alternate form 

v £ (t) = e “ i£Qt w £ (t) (39) 

where 

2 TT 

W £ (t) = — | wu,t) e ix<j) d<j) . (40) 
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The w £ (t) deflections, defined at the reference blade posi- 

ion as it rotates, will be the dependent swashplate unknowns from 
this point on. Once these quantities are found, the actual swash- 
plate displacements can be found from either of the following two 
equations : 

rotating frame 

°o ^ 

w(<i>,t) = l e l£< ^ w £ (t), (41) 

£=- “ 

non-rotating frame 

v ( 0 , t) = I e i£(0 " ftt) w (t). (42) 

l = -00 


Substitution of the expression for (t) as given by Equation 

(39) into Equation (37) and multiplying by the factor e l£fit yields 
the governing equation for each w^ (t) which is 


3 2 w a ( t) 

at 2 


3w. (t) 

2i4-fi — £ 2 £2 2 w, (t) 

a t: & 


+ SiH. w Cti 

_ Q X/ 

R d 
sp 


= e i£f2t [q m - a^t). + S, 3 £ (t) 

sp 


(43) 


Definition of applied loading function .- The harmonics of 
swashplate loadings required by the right hand side of this equa- 
tion can be obtained in terms of the swashplate displacement, 
twist, bending slope and collective base deflection. 


Recall that Q. (t) was obtained from the Fourier transform of 


Equation (23) so that 


.2 IT 


(t) = 27 


Q(0,t) e 12,0 de 


(44) 


and similarily from Equations (24) and (25) 

r 2ir 

i 


= 2 TT 


a (0 , t) e ^ £0 d0 


(45) 
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and 


V t! 


i 

2t r 


2TT 


6(6, t)e 


-l£0 


de 


(46) 


One contribution to Q (0 ,t) comes from the discrete forces due 
to the rotating control rods which couple the swashplate with the 
rotor blades. The other contribution comes from the forces in the 
spring-damper units shown as item (2) in Figure 6. 

A Dirac delta representation of the concentrated forces, 
moments, and torques leads to the following expressions for 
6 ( 6 , t) , Q ( 0 , t) and o (0 , t) ; 


N,. 


3(0, t = l 


m=? 


p (t) as (e 

m m 


H't $ +(().) 

m Y mp 


N 


e. s. 


+ 1 6(6 ~ X-)b. <jk^ + ~ ) a(0/t) 


j=i 


3 3 


N 


- i 6(0 - X.) (k i + c Itt) 4 .( 0 , t ) 

j-1 3 *-j 9t 


N, 


Q ( 6 , t ) = l P ( t ) 6 ( 0 - fit - <j> + 4) 

m TUI TV 


m=l 


m mp 


N 

e.s. 

I 

j-1 


X ) (k. + c -r) d(0,t) 

3 J ] ^ 


(47) 


(48) 


and 


o (0 , t) 


N 

6 • S • 

- j 6(0- Xj ) (k Q> + c ft) ®(e.t) 

j-1 3 3 


( 49 ) 


where d(0,t) 


v (0 ,t) 


b^$ (0 , t) 


u 0 (t) , 
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is the number of blades, 

P (t) is the force in the m— control rod, 

d m is the rigid offset of the m— — control rod attachment 

point from the neutral axis of the ring, positive outward, 


<f>m = — - for equal spacing between control rods, 

N b 

4> m p = blade wrap-up angle, 

c is the number of spring-damper units directly 
supporting the swashplate, 

"ttl 

X . is the azimuthal angle locating the j — support 
3 relative to swashplate x,- axis 

til t 

b . is the offset of the j— — linear spring attachment point 
from the neutral axis of the ring, positive inward 


k. and c. are the linear stiffness and damping values for 
3 3 

the j— spring-damper unit, 


k Q and c Q are the torsional stiffness and damping values 

U . D . 

3 3 

the j' — torsional spring-damper unit counteracting local 
lateral swashplate rotation, 


k, and c, are the torsional stiffness and damping values 

♦j b 

for- the j — torsional spring-damper unit counteracting 
local longitudinal swashplate rotation. 


u 0 (t).is the vertical downward displacement of the base 
which supports the spring-damper units, 

and 6 ( © — 9 o ) is defined such that 


•Zn 

f (6)6(8 - 0 o ) R d0 = f ( © o ) (50) 

J sp 

0 

(0 O being an arbitrary azimuthal angle). 
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The variable, <j> , is further defined as the azimuthal location 
in th 

(phase angle) of the m-^— blade relative to the swashplate x^-axis 

when t equals zero and $ is taken as zero in value. 

^ mp 


The required expressions for 0^ (t) , Q^(t), and a (t) , needed 

in Equation (43), are obtained by substituting the forms for 0 ( 0 , t) 
, Q ( G , t) , and a ( 0 , t) , as given by Equations (47) , (48), and (49), 
into the integrand of Equations (46) , (44) , and (45) , respective- 
ly. Taking proper care with regard to the integral of the Dirac 
delta functions, the Fourier transforms of 0 ( 6 , t) , Q ( 0 , t) , and 
and a ( 0 , t) are 


N. 


B £ ( t ) 


1 — Y P (t) d 

P ^ _ -m v 1 n 
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2 7T R S m ' ‘ "m 
sp ra=l 


N 


e. s . 


e. s . 


1 w • ^ * 1 

— ,L % + c * j 3t> e 


Q,(t) 


°t (t) = 


2 7T 

R 

sp 


] 

1 

2 7T 
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■ R .4 

sp m= 
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2 TT 

R sp 


1 

2 7T 

R sp 


— i£x * 

:) bj d (xj ft) e 3 


N 


-iMat + + m - v p > 


‘ m 


N 


e,s. 

r» 

i 


(kj + O. 3 £) d( Xj ,t> e 


-i*X. 


N 


e . s . 


I (k 6j + c ej ft* ^Xj'** e 




(51) 


(52) 


(53) 


where d( x .,t) = v( x .,t) - b. *( X .,t) - u n (t) . 


The quantities; v( X j,t), 0( X j,t), $( x ^,t) and d( x ^,t); in the 

previous equations . represent the swashplate? local deflection, 
local bending slope, local twist and local spring-damper unit 

deflection at the azimuthal location, x ^ , of the j — spring- 

damper unit. These variables can be expressed in terms of the 
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Fourier harmonic components of the swashplate deflection in the 
rotating frame. 


Replacing 0 with the specific value, Xj and altering the 

dummy variable over which the summation is taken. Equation (42) 
becomes 

°° ip(x . -fit) 

v(x^t) = l w (t) e 3 (54) 

J p=-00 

Using the same procedure with the form of Equation (29) , 

0 ( X j f t) can be expressed as 


00 IPX - 

© ( X • r t) = l 0 ( t) e 3 (55) 

J p— — CO ^ 

Replacing SL with p in Equations (39) and (37e) , the equations 

v (t) = e _ipfit w p (t) (56) 

and 


0 (t) = ip v (t) / R 
P F P sp 


are obtained. Substitution of Equation (56) into the last equa- 
tion results in an expression for 0 (t) which when inserted into 

Equation (55) yields the expression 

» ip(X4" t) 

0 (Xj i t) = « — I ip w_ (t) e 3 (57) 

3 sp p--°° v 



The quantity $(xj/t) 
,t) in Equation (55) , 


can be expressed in the same form as 
that is 


00 ipX . 

$ (X.,t) == l $ (t) e 3 

J r^=r— »°° t' 


(58) 


Altering the dummy variable, l, to p in Equation (37d) and drop- 
ping the 8 (t) term since its contribution, based on a comparison 

of the effects of the terms on 5> (t) , is negligible results in 

Jr 

* P (t) = s p v p (t) 
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Substitution of Equation (56) into this equation yields an ex- 
pression which when inserted into Equation (58) yields the ex- 
pression for $(Xj/t) in the form 

< Kx-/t) = l S w (t) (59) 

“r)=— -oo ¥ ir 


By substituting Equations (54) and (59) into the expression 
for d(Xjft) following Equation (53), an expression for d(Xj,t) as 

a function of swashplate rotating frame deflections and the collec- 
tive base plate deflection can be obtained in the form 




ip (x.-fit) ' 

(1 - b.S ) w (t) e 3 
DPP 


" u o 


(60) 


The collective base deflection, u Q (t) , can be related to the 

deflection of the swashplate by consideration of the force equili- 
brium of the collective base plate which results in the expression 


(K + C |^) u o (t) 


l (k. + c. |x) d(x v t) 

j=l : 3 dt 3 


(61) 


Insertion of Equation (60) into this expression provides a relat- 
ionship between the base plate and swashplate deflections in the 
form 


N 


e . s , 


® + c It) + T' < k i + c i ft> 


U 0 (t) 


N 

°° e s * 


%L. ih u - b jV ft 3 * C J # w P (t) e 


ip (x .-fit) 


(62) 


which cannot be solved for u 0 (t) without use of Laplace trans- 
forms , 


Expressions relating the swashplate loading harmonics to the 
swashplate deflections in the rotating reference frame and the 
collective base deflection can now be obtained. Insertion of 
Equations (57) , (59) , and (60) , as required, into Equations (51) 
through (53) results in the equations 
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3 3 


4 , i, ^ 0 1 1 0 at' o 

p=_co j = l j 0 j a 2 IT R gp 2 


(65) 


Application of Laplace transforms .- The real-time swash- 
plate governing equation is Laplace transformed to replace the 
derivatives involved with algebraic forms in addition to pro- 
viding compatibility with the blade equations which will be de- 
veloped. In addition to the specialized Laplace transformations 
which were presented previously pertaining to the development of 
analysis for representation of the gyroscope control system, two 
more are required for this section. These are 

(a) L^ (t) e = Wp (s-i£fi+ipO) 
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(b) L t ( w ^ (t) e - ( s -i£0) w^ (s-iiLfl+ipfi) . 

The Laplace transformed version of Equation (43) assuming 
quiescent initial conditions is given by 


y (s? 2 - 2i2,fis - Z 2 Q 2 \ + £liiL w . (s) 

p 3 £ 


T riiflt ^ . v ^i&fit ii „ , ,, , „i£nt „ 0 ,.r\ , cc , 

L t e Q Z ■ ~ e R — £ ^ + e S £ • (66) 


The Laplace transform of each of the terms appearing on the right 
hand side of the above equation can be obtained. 

Multiplication of Equation (63) by e 1Jifit and application of 
the Laplace transformation yields the expression 


N. 

L t[ ai<!!!t V«] - | 
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J p l 2tt R j 

■ sp 


e*s„ u 0 (s-i£Q) 

+ X Kc <i TTr e ' 

3=1 sp 


where for convenience 


KC = k, 4- (s-iaft") c . 
a 3 3 


It is important to realize that the derivatives involved in Equa- 
tion (63) through (65) do not operate on the multiplying factor, 

e x£ftt ^ rph e quantity, w^ (s-i£fi+ipfi) , is not to be interpreted as 

a product but rather as a shifted argument of the function where- 
as the term, (s-i£fi) , does represent a product of two terms. 


In a similar manner as was used to obtain Equation (67) , 
Equation (64) and (65) can be used to obtain the expressions 
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where for convenience 


and 


KC, = k, + (s-iafi) c , 

■9 9 • 9 


KC 0 = k 0 + (s-iafi) c Q 
a j j 


The Laplace transformed version of Equation (62) makes it 
convenient to solve for uq (s) as 


u 0 (s) 
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“e.s. ipx* _ 

I I (1 - b.S ) KC 0 e J w (s+ipO) 
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The expression for u Q (s-i£Q) which is required in Equations (67) 
and (68) may be obtained by replacing s by (s-ijifi) ( shifting the 
argument of the function)" in Equation (70) which results in 


u 0 (s-iJlfi) 


N 

” e.s- IPX-; _ 

l l (1 - b.S) KC p e J w^ (s-ijiQ+ipo) 

p=-co j = i 0 P 
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e.s. _ 
K + (s-itfi) C + l KG 

j=l 


This expression can be substituted into Equations (67) and 
(68) which with Equation (69) can replace the terms appearing on 
the right hand side of Equation (66) resulting in the final 
swashplate control system governing equation of motion. By re- 
placing the index l with q and then the index p with t in the 
resulting expression, for convenience in later use, and multi- 
plying by 2 7 tR , the following system of coupled equations for the 
sp 

Laplace transformed versions of the swashplate unknowns, w^'s, in 
terms of the control rod forces, .P (s') , results. 
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where in is the mass of the swashplate and is equal to 2 ttR y, 

sp 


When the governing matrix equations for the rotor blades are 
developed in the next section it will be the control rod forces 
(P^Cs)) which couple the swashplate equations or the gyroscope 

equations to the blade equations. These relations will be used to 
develop the final complete set of equations for the control system- 
coupled rotor blade system. 


Representation of Rotor Blades 
by Transfer Matrices 

In this section the analysis, using a transfer matrix 
approach, is developed for representation of the blade behavior 
due to forces and moments of aerodynamic, mass, and inertia origin 
acting on the blade. Since the transfer matrix method is compli- 
cated by consideration of aerodynamic blade sections , the matrix 
approach is presented, initially, for sections not involving aero- 
dynamic considerations and then extended to the form required for 
aerodynamic representation. The resultant analysis provides the 
necessary representation for blade behavior including coupling 
with the control system behavior. 

General blade model. - A lumped parameter approach has been 
utilized as the means of representing the allowed blade distributed 
properties which were listed previously in the section pertaining 
to helicopter rotor system modeling. Each rotor blade is repre- 
sented by a specified number of basic sections. Each blade 
section can have either (1) concentrated mass and inertia, uniform 
elastic properties, and geometric characteristics such as offsets, 
precone, presweep, and blade twist, or (2) aerodynamic and certain 
geometric characteristics. The separation of aerodynamic charac- 
teristics, section type (2), from the characteristics of section 
type (1) is necessary for proper determination of the modified 
aerodynamic forces acting on the blade section when the free-wake 
coupled solution method is to be applied. However, if the free- 
wake coupled iteration solution is not to be applied to the model 
of interest, the model can consist of sections each having both 
type (1) and type (2) characteristics. The general construction 
of a blade section showing the order in which specific character- 
istics are considered in going outboard to inboard along the 
section is illustrated in Figure 8. 
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concentrated 
mass and inertia 


elas tic 
section 



elastic axis 
outboard 


Figure 8. A general blade section 

It should be noted that concentrated spring-damper units , bends , 
and rigid offsets are allowed in three mutually perpendicular 
directions and that the section center of mass and elastic axis 
may be offset from the quarter chord (pitch axis) of the blade. 

Coordinate system involved. - The coordinate systems defined 
for representation of blade behavior consist of a fixed (non- 
rotating) coordinate system, rotating coordinate systems, and local 
blade coordinate systems. The orientation of the axis of the 
fixed shaft (hub) coordinate system has been specified previously 
in the discussion concerning the swashplate coordinate systems. 

The origin of this system is taken to be on the shaft centerline 
in the plane of the rotor blades when the blades are not preconed 
or preswept. The shaft (blade) rotating coordinate systems have 
their x and y axis in the same plane as the respective axes of the 
non-rotating shaft coordinate system but are rotating at a speed 
of fi relative to the fixed shaft system. For each blade there 
exists a rotating shaft coordinate system whose x-axis is located 
along the spanwise axis of the blade when in the unconed and un- 
swept position and which can be related to the shaft fixed co- 
ordinate system by the coordinate transformation. 
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1 rm 
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f > 
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(72a) 

^rrrr 


0 

0 
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where fit and <J> have been defined. 

Y m 

The blade local coordinate system C x rb ry rb / z rb ) is obtained 
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rm 


Figure 9. Blade local coordinate system and state 

variable orientation, slash marks refer 

to x . , y , , and z . blade coordinates , 
rb -^rb' rb 

respectively 



from tile associated shaft rotating coordinate system 
^rm ,y rm ,z rm^ by performing three orthogonal rotations; $, 6, and 
in consecutive order where 


$ corresponds to forward presweep 

0 corresponds to downward precone, and 

y corresponds to nose-up pre-twist and/or 

collective angle of attack 


The relationship between the blade local coordinate system and its 
associated shaft rotating coordinate system is depicted in Figure 
9. The mathematical relationship between the two coordinate sys- 
tems is represented by the coordinate transformation 
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where a short form of denoting the sine and cosine functions 
has been utilized, Ge.g., .S$ represents sin$ and C$ represents 
cos$) . Since the transfer matrix approach which will be de- 
veloped in the next section relates the blade deflections, 
slopes, forces, and moments along the blade relative to its 
local coordinate system, the transformation shown above can be 
used to obtain blade variables in the disk plane. 

Transfer matrix representation of blade model without aero- 
dynamic loading and hinge or bearing discontinuities.' - The trans- 
fer matrix approach utilized for a blade without aerodyn ami c s 
included in the blade model is presented initially to provide a 
mathematical foundation for the presentation of the modified 
transfer matrix approach required for representation of aerody- 
namic effects on the blade behavior. The basic transfer matrices 
to be utilized in this section consist of essentially three basic 
types. These are (1) section individual transfer matrices which 
are used to construct- (2) section transfer matrices which in turn 
are used to construct (3) associated transfer matrices. Corre- 
sponding to these types of transfer matrices there also will be 
involved corresponding types of forcing functions. Each type of 
transfer matrix and forcing function can be represented in a real- 
time form which upon application of Laplace transforms results in 
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a related transfer matrix involving complex variable notations. 
In this initial analytical development for representation of the 
rotor blade behavior it is convenient to consider the basic 
real-time relationships only briefly and carry out the develop- 
ment of analysis utilizing Laplace transformed representation. 


A section transfer matrix can be developed which relates a 
state vector at either end of a section such that 


(D- 

m 


= . 6 ]» i) H» il+ 


l(i) 

m 


(73) 


The + and - signs in Equation (73) refer, respectively, to the 
outboard and inboard ends of the ith section of m th blade, the 
sections being numbered outboard to inboard. The state vector. 


S>, is a column of twelve quantities expressed as 


; S ) = <u x N T u y M z 


V - u 4> M V 

y z y y z 


where u y , u z 


N f V y? V z 


V *y' 


are blade displacements in the blade local 
x, y, and z directions; 

are a^ial force, edgewise shear force, and 
flapwise shear force, respectively; 

are local torsional deflection, flapwise 
bending slope and edgewise bending slope, 
respectively; and 


T, M , M are local pitching moment, flapwise bending 
y z moment, and edgewise bending moment, 
respectively. 


The positive sense of the twelve state variables comprising the 
state vector has been shown in Figure 9 as they would act on the 
outboard end of the next inboard section from the section to which 
the transfer matrix has been applied. The vector, (F), in Equation 
(73) is a section individual forcing function consisting of mass, 
inertia, and maneuver loading effects on the shear and moment 
variables in the state vector and may also be considered to be the 
section forcing function since aerodynamic considerations are 
being omitted initially* Contributions to the section forcing 
function due to other individual characteristic portions of a 
section do not occur. 


In general the transfer matrix is a 12 x 12 array whose ele- 
ments are differential operators when each of the state variables 
is a real time quantity. By taking the Laplace transform of 
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Equation (72) the differential operators are converted to complex 
linear algebraic equations. Assuming quiescent initial conditions 
the transformed version of Equation (73) is 



The section transfer matrix 



can be expressed in terms 


of more fundamental transfer matrices which are associated with the 
mass, elastic, and geometric properties of the section, individual- 
ly . Thus 



where the individual transfer matrices are constructed such that 

H accounts for a rigid offset in the elastic axis or a 
translation of the local coordinate system; 

H accounts for a bend in the elastic axis or more generally 
a rotation of local coordinates; 

E accounts for a uniform elastic section with centrifugal 
L--! stiffness included; 

SK accounts for a concentrated spring-damper (used normally 
*- -I for lead-lag hinge only) ; and 

[a] accounts for a concentrated mass and inertia. 


Taking the Laplace transform, it can be shown that 



where each matrix consists of time independent terms and/or Laplace 
operators. In practice some of the individual transfer matrices 
may not be necessary in representing a blade section and would be 
omitted from Equation (76) . 


In order to develop the associated matrix across a rotor blade, 
consider the transfer across the (i+l)th blade section which from 
Equation (74) is given by 


sl ” 

m 


( i+i ) isi ^ i+i ) + + 

m m 


59 



However , 



(i+1) + 



(78) 


Substitution of Equation (74) into Equation (78) and insertion 
of the resulting expression into Equation (77). yields, the 
equation 



(i+1) - 
m 


yi (i+i) 

L> 



+ 



(i+1) 

m 



+ 


Mm 1+1> ‘ 79 > 


for transferring across the (i+1) th and ith blade section. Extend- 
ing Equation (79) , the expression which reTates the state vector 
at the inboard position of the ith section to the state vector at 
the blade tip can be written in tEe form 




i-1 

+ I 

j=l 


[k=j+l 


fsl « |f} { ^ 

_Jm ( Jm 



The blade tip state vector denoted by 



has the 


(80) 


moments and shears equal to zero due to the free boundary condition. 
It is convenient to use this fact in the blade transfer process 
and define 
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This matrix equation can be rewritten as 


rl (Tip) 

7 m 



( 82 ) 


The reason for writing Equation (81) in the form shown is connected 
with the modified transfer matrix procedure used to solve the final 
system of equations for the dynamic forced response. 


The fact that the 


VI Ti P 

Jm 


matrix is a 12 x 6 instead of a 


12 x 12 reduces the number of computations required to obtain the 
associated transfer matrix for the whole blade. 


On insertion of Equation (82) into Equation (80) , the 
equation 






+ 



(83) 


is obtained where the associated matrix is defined as 




(84) 


and the associated forcing function column is defined as 




j=l 


i 

n 

k=j+l 



CD 

m 


(85) 


Equation (83) is the blade matrix equation which is valid for any 
section which is outboard of all possible hinge and bearing dis- 
continuities in the blade model. 

Representation of blade hinge and bearing discontinuities 
without aerodynamics. - The analysis being developed includes the 
representation necessary for consideration of flap hinge, pitch 
bearing, and rocker arm attachment point discontinuities in the 
blade model. At these locations there are discontinuities in flap 
angle, torsional deflection, and blade torque, respectively. 

The state vector relationship across the pitch bearing dis- 
continuity location may be written in the form 
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- W£ ea (s6 > 

where the column vector 

0 0 0 0 0 \ . ( 87 ) 

J m 

The state vector relationship across the flap hinge discon- 
tinuity location, provided there exists a flap hinge in the blade 
model, is given by 

Mr apl 

where the column vector 

{l}^ lap = jo 0 0 0 0 0 0 0 0 IJ y 0 0 } m . (89) 

In transferring across the rocker arm attachment point it is 
assumed that only the torque applied to the blade from the control 
system is significant. Thus the transfer inboard across the cont- 
rol torque point can be written as 


- 


s}' £lap,+ 


flap 

m 


(88) 




(c. to- 


rn 


, (c.t.) + _ Mc.t. 
r m r m 


(90) 


where the column vector 


m 


000 AT 00000000 


m 


(91) 


The quantities AT, A , and are additional blade unknowns which 

are to be solved for in the final system of equations. 

Equations (90) , (86) , and (88) can then be replaced by the 
matrix equations 


sl * c * t * * ” 

m 


s Hc.t.) + 
in 


m Jm 




(92) 


( 93 ) 
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{gjUiap)- . |s}« la P> + - {i}f ap {x}* flap 
where the following arrays only have one element: 



( 94 ) 


(95) 


(96) 


Expression (83) can now be replaced by the matrix equation 
which provides the state vector at the inboard end of an arbitrary 
blade section which is 



where by letting N c t , N^ ea , and N f^ a p be the blade section numbers 

immediately inboard of which the control torque, pitch bearings, 
and flap hinge are respectively located, the new matrices appearing 
in Equation (97) are defined as 




i < 


N 

C. 


t. 


(98) 
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i > N 


fea 


rlfea 


4 h- 


1 = N 
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1 < Nfea (99) 



A 



1 < N f lap (100) 


With the above definitions the state vector, jsjj^ ”, when evaluated 

at i equal to N c t , N^ ea , or N f^ a p will be immediately inboard 

of the control torque point, the pitch bearings or the flap hinge, 
respectively. 


Application of blade root and discontinuity boundary condi- 
tions when no aerodynamxc loading. — By use of the boundary 
conditions which occur at the blade root and the conditions which 
must occur at the control torque point, pitch bearing, and flap 
hinge, the governing blade equations in terms of blade tip un- 
knowns and discontinuity unknowns can be obtained when aerody- 
namic sections are not allowed. 


Displacement and slope boundary condition applications: 
Letting NS be the number of blade sections, the state vector at 
the blade root is determined from Equation (97) by letting i 
equal NS. If the boundary conditions at the root of each blade 
are such that all displacements and slopes vanish, then these 
conditions are expressed by the matrix equation 



U0.D 
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Or, from Equation (.97) 


m 




* (Tip) 

m 




where 


- >]{ a }™ S! H; flap + HHr 1 = w 



1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 


(103) 


Equation (102) provides six equations for each blade. However, 
there are nine possible unknowns due to the control torque, pitch 
bearing, and flap discontinuities. The other three equations are 
obtained from the condition of zero bending moment and torque at 
the flap hinges and pitch bearings respectively, in addition to the 
relationship between the control torque and the blade torsional 
deflection at the rocker arm attachment point. 


Application of zero torque condition at pitch bearing: 
Consider the condition of zero torque at the blade pitch bear- 
ings. The state vector at the pitch bearings is given by 

M (N fea } 

{S> m . In order for the torque (which is the fourth 

quantity in the state vector) to vanish the equation 

r\ (N f _ 

( r -){ s }m = 0 (104) 

must be satisfied. Or equivalently 






= 0 (105) 



where 


( r 


») 


(0 0010000000 0 ). 


Application of control torgue condition for coupling with 
swashplate or gyroscope control system: The equation for the con- 

trol torque must reflect the coupling between the blade and the 
swashplate or gyroscope motion. In real time the balance of 
forces on the mth control rod which connects the rocker arm to the 
swashplate is given by 


(AT) = 
m 


a k (1 + t tt) 
m m m 8t 


(N . ) ~ -<t> ) 

a U ) c,t * - l (1 + d S.)w 0 (t) e m mp 
m Y x m . L m 2 . 5 

2 .=-» 


(106) 


where k and t are the linear stiffness and damping coefficients 
m m 

of the mth control rod, respectively, and a m is the distance of con- 
trol rod rocker arm attachment point aft of quarter chord. Divid- 
ing by k m and taking the Laplace transform yields 
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(107) 


Referring to Equations (90) , (91) , and (92) 
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and also 
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where |^r 3 j = (0 0 1 0 0 0 0 0 0 0 0 0 ) . 

Thus, the control torque equation can be written in terms of the 
blade tip and internal discontinuity unknowns as 
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( 110 ) 


If the swashplate is replaced by a gyroscopic control system, 
the equation for the control torque must reflect the coupling 
between the rotor blade and the rigid body motion of the gyroscope, 
In real time the balance of forces on the mth control rod which 
connects the rocker arm to the gyroscope is given by 


(AT) = a k (1 + t .V) 
in in xn m at 


(w , ) 

a (4) ) c,t - 
m Y x m 


I w'(t) 
£=-l ^ 


(HD 


where the w^ quantities now stand for the rigid body motions of the 

gyroscope (i.e. , see Equations (17), (18), and (19) for equations 

of motion of the gyroscope) and are given as 
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( 112 ) 


Dividing by k m and taking the Laplace transform yields 
i (N ) 1 

kl CATV - a 2 (1 + IS) (?) c * t * 


m 


m m 


m x m 


+ l a (1 + t s) wf (s) = 0 . 
£=-1 m m * 


(113) 


Therefore, whether the rotor blades are coupled through the swash- 

plate or the gyroscope, the general formulation is essentially the 

same with the only difference being the functionality involving the 

(<!> - <f> ) term, 

m mp 


Application of flapwise moment condition at flap hinge: If 

a flap hinge is present the condition that the flapwise moment 
vanish gives the last equation needed. Since the flapwise 
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moment is the eleventh quantity in the state vector, this condi- 
tion is satisfied provided 


!N flap> 


where 


( rii )H 

jr n J = ( 000000000 0 10) . 


( 114 ) 


Referring back to Equation (97) , the zero flap hinge moment is 
expressed by 


bo[C flapi {s}: <Tip) - 


( Wj* 

m 




flap 
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(. 115 ) 


For a blade without aerodynamic loading, Equation (97) re- 
lates the state vector variables at the inboard end of the i th 
section of the mth blade to the blade tip unknowns and discon- 
tinuity unknowns and thus may be utilized to obtain the blade 
shape if all unknowns have been determined. Values for the un- 
knowns are obtained by use of Equations (102) , (105) , (115) , and 

either (110) or (113) depending on the type of control system 
involved. 


Shift of transform variable in sectional transfer matrix 
Equation and its interpretationT - Prior to the inclusion of aero- 
dynamic blade sections into the blade model which compounds the 
complexity of the transfer matrix approach, consideration of the 
representation implied by the use of Laplace transform shifted and 
unshifted variables will be advantageous. 


In the previous development of the transfer matrix equation. 
Equation (74), the Laplace transformed arrays are a function of the 
Laplace transform variable, s, although this dependency was not 
explicitly shown. This dependency is directly observable by re- 
writing Equation (74) in the form 

{® ts) }m i)_ ‘ [ 5(s) ]m il { I(S) }m i)+ + { F(s) ' (116 > 

In the dynamic response type of solution the total state vectors 
and forcing function corresponding to those in the above express- 
ion are expressible in real-time notation in a Fourier series form 
exemplified for the state vector by 

jstt)}' 1 *-- l { S "If - 
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where the S n harmonic coefficients may be complex variables and 
correspond to the commonly denoted; 0/rev. , ± 1/rev. , and etc.; 
blade responses and forcing functions for the state vectors and 
forcing function arrays, respectively. Substitution of iNpg for 
the Laplace transform variable, s, in Equation (116) provides the 
transfer matrix equation relating Np/rev. coefficients of the state 
vectors on both ends of a section and of the forcing function on 
the section. 

This can be shown by first applying the Laplace transform 
operations as specified in Equation (116) to the real-time state 
vectors and forcing function such that 


f s n 1 (i) 

\ (s-i'nfi) / m 



f S n 1 ( i ) + 
\ (s-inn.) J m 
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results. Multiplication of both sides of the above expression by 
(s-iNpft) yields the equation 
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n =-co 1 (s-inf2) jm 

n/Np 


which in the limit as s approaches iNpn reduces to the form 



relating the Np/rev. harmonic coefficients of the state vectors 
and forcing function involved.. Thus, substitution of iNpfi for s 
in Equation (116) corresponds to the relationship shown in the 
last expression. 


Shifted Laplace transform variables of the form (s-ikfi) re- 
present harmonics relative to the main harmonic of interest de- 
noted by the value of s. By shifting the argument of Equation 
(116), that is, replacing s by (s-ikO) , the expression 

js (s-lkfl) = jG(s-ikfi)J^ l) |s(s-ikn) |^ l)+ + |f (s-ikQ) (117) 
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is obtained. Application of the same procedure that was used on 
Equation (116) to Equation (117) yields when s is replaced by iNpfi 

r s N P -k} ai - = j- g(1(Hp . k)a) -|a)| s Np-k|a) + + p P -k|«i, 

which relates the (Np-k) /rev. harmonic coefficients of the state 
vectors and the forcing function. In this equation and Equation 
(117) a positive value for k corresponds to a downward shift in 
the harmonic from that specified by s whereas a negative k de- 
notes an upward shift in harmonic. The quantity, s, taken as iNpO 
in Equation (117) provides the relationship shown in the last ex- 
pression. This shifted Laplace transform concept is valid for 
any originally time dependent variable used in this report. 

For convenience in the representation of the transfer matrix 
approach when aerodynamic sections are included in the blade model 
Equation (117) can be written in the subscripted form 



(118) 


Transfer matrix representation of transfer across an aero- 
dynamic section .- The transfer across an aerodynamic section in 
real-time notation can be represented by the matrix relationship 


(i) = c f s \(i) + + D (i) 

J dt ( J _ 


)U<i) + 


+ id 


(i) 


(119) 


where H ^ and H ^ are transfer arrays involving periodic 

' (i) 


functions and the column vector, |d j / is the aerodynamic forc- 

ing function also consisting of periodic functions. Each of these 
arrays can be expanded in a Fourier series form. For example, 


H (1) " JJJ " 


) inflt 
e 


(120) 


Defining jXJ ^ ^ relative to the initial position of the first 
blade and replacing Qt with (Qt+^ m -^ m p) , Equation (120) becomes 


for the mth blade 


tr ■ i-fei- 1 


\ in (Ot+<j> — 4> ) 

) m mp' 


( 121 ) 


Applying the same procedure to the other two arrays appearing in 
Equation (119) and substituting into Equation (119) results in 
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The Laplace transform of Equation (122) , assuming quiescent 
initial conditions and also shifting the Laplace transform 
variable, s, is 
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(123) 


The summation on the first term of the above expression can 
be taken from n=-N f to where is the number of harmonics 

above and coupled to the frequency being solved for. The summa- 
tion on the second term can be taken from n=-N^^ to where 

is maximum number of aerodynamic forcing function variable har- 
monics not including the steady, that are obtainable. 

However, when these limits are placed on the respective 
summations involved in Equation (123) , it can be observed that the 
equation representing shifted harmonics of the state vector, k 
not zero, involves harmonics further removed from the main harmonic 
of interest. For example, if N f is taken as 1 corresponding to 

inclusion of interharmonic coupling 1/rev. above and below the 
harmonic corresponding to k equal to zero, the version of Equation 
(123) for k equal 1 would include a term involving (n+k) equal to 
2 which corresponds to a harmonic 2/rev. below the main harmonic 
and similarily, for k equal -1, a term involving (n+k) equal to -2 
results. This addition of extra harmonics above and below the 
range specified by the value for N f is not desired since these 

additional harmonics would in turn have to be represented and 
thereby result in further additional harmonics and so on. 

To overcome this problem the summation involved in the first 
term of Equation (123) can be taken from n=-Nj-k to N f -k. How- 
ever, since n includes (-k) the same result is obtained if the 
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summation is taken over n=-N f to N f and n in Equation (123) in 
the first term is replaced by (n-k) . This yields 


kjm 
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which can be written in the form 
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(125) 


where Q'n-k represents the aerodynamic section transfer 

matrix for a given section with (n-k) denoting which harmonic of 
the C and D matrices is involved and n representing the shift on 
s inherent in the matrix. 


Transfer matrix representation of blade model with aerody - 
n ami cs but without hinge or bearing discontinuities.- Equation 
(118) representing the transfer across a non- aerodynamic section 
and Equation (125) representing the transfer across an aerody- 
namic section provide the fundamental relationships necessary to 
construct the associated transfer matrix representation relating 
the state vector at the inboard end of a blade section to the 
blade tip unknowns. 

In order to include the aerodynamic interharmonic coupling 
effects, Equation (83) is modified to the form 
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where the k and n subscripts denote the degree of harmonic shift 
relative to the harmonic specified by the Laplace transform 
variable, s, as previously noted. The modified associated 

transfer matrices of the form, , represent the effects of 


n-shifted harmonic coefficients of the blade tip slopes and de- 
flections on the k-shifted harmonic coefficients of the state 
vector variables at the inboard end of the ith section and for 
the mth blade. The modified associated forcing function vectors 


of the form, 



(i) 

m 


represents the effects of the k-shifted 


forcing function harmonic coefficients on the k-shifted harmonic 
coefficients of the state vector variables at the inboard end_ 
of the ith section due to aerodynamic, mass, and inertia loadings 
on and outboard of this blade section. 


Since inclusion of aerodynamics into the blade model results 
in interharmonic coupling which is apparent in Equation (125) , 
the subscripted associated transfer matrix arrays and associated 
forcing function vectors cannot be expressed in the simple forms 
exemplified by Equations (84) and (85) when aerodynamics are in- 
cluded. Instead, the relationships between the associated 
transfer arrays and forcing functions at the inboard end of a 
section and the associated transfer arrays and forcing functions 
at the outboard end of a section for the two types of blade 
sections can be determined. 


For a non- aerodynamic section. Equation C126) with the (i) 
superscripts replaced by Ci-1) superscripts can be inserted into 
Equation (118) for the (i)+ superscripted state vector (making use 
of the form of Equation (78) ) such that the equation 



results. From a comparison of this expression and Equation (126) , 
the relationships of the associated transfer matrices and asso- 
ciated forcing function vectors of a section to those of the 
previous section can be expressed as 



(127) 

(128) 
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For an aerodynamic section the same form of comparison can 
be used. In this case, the superscript in Equation (126) is 
altered in the previous manner but the k subscript is replaced 
by p such that when the resulting expression is inserted into a 
form of Equation (125) where the dummy index n has been replaced 
by p in the first term the following equation results 
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Referring back to Equation (126) 
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It should be noted that the associated transfer matrices 
for k not equal to n consist of all zero elements until an 
aerodynamic section is encountered since non— aerodynamic stations 
do not provide interharmonic coupling. This can be observed by 
noting that the definition of tip state vectors in terms of the 
tip unknowns by application of tip boundary conditions is repre- 
sented by a subscript form of Equation (82) which is 
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, (Tip) 
m 



(Tip) 


where the matrix 


W Tip i 

L k J m 


is defined for each value of k as pre- 


viously defined by Equation (81) . The above expression can be 
extended to the form 



(132) 


where the 



(Tip) 

m 


•array is defined 


for k equal to n as in 


Equation C81) and consists of all zero elements if k is not equal 
to n. Thus, if the first section is a non-aerodynamic section. 
Equation (127) would be replaced by 



based on a subscript form of Equation (84) when i is equal to 1. 
Thus, from this expression and Equation (127), the associated 
transfer matrices for k not equal to n can be seen to consist of 
all zeroes until an aerodynamic section is crossed at which 
point the associated transfer matrices for k not equal to n re- 
ceive a contribution due to interharmonic coupling as shown by 
Equation (130) . 

With the inclusion of aerodynamic sections. Equation (126) 
has become the general matrix transfer equation for sections 
outboard of hinge or bearing discontinuities where the asso- 
ciated transfer matrices are related by either Equation (127) or 
(130) and the associated forcing function vectors are related by 
Equations (128) or (131) 

Hinge and bearing discontinuities representation modifica- 
tions due to inclusion of aerodynamics.- The state vector at the 
inboard end of any arbitrary blade section on extending Equation 
(97) to a subscripted form is defined by 
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(133) 


Since, normally, an aerodynamic section would not occur inboard of 
the flap hinge, pitch bearings, or rocker arm attachment point 
Equations (98) , (99) , and (100) can be extended to 


{' 


k,n 
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< N fea + l> 
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(134) 


where 


( X k,nJ 


•t. h ] 

' rk,nj 


f ea 
m 


and 




flap 

m 


are defined in 


Equation (96) when k is equal to n and consist of all zeroes if k 
is not equal to n. If an aerodynamic section does occur inboard 
of the flap hinge, pitch bearings or rocker arm attachment point, 
then, as the result of a transfer across this section, the arrays 
defined in Equation (134) are altered such that, for example, 


k,njm 


(i) 


pI-N (rtl) {kp' n l m 


(i-i) i( p- k) ( w 


(135) 


and similar equations occur for |c^ and |d^ n 


(i) 

m 


Application of blade root and discontinuity boundary condi- 
tions with aerodynamics. - From the general form for the state 
vector inboard of any arbitrary section as represented by Equa- 
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tion (133) the boundary conditions at the blade root and the dis- 
continuity conditions can be applied as in Equations (101) 
through (115) to yield the representative equations when aero- 
dynamic sections are included. 

Displacement and slope boundary condition applications: With 

the slopes and displacements being zero in value at the blade 
root irrespective of- the value of k, the insertion of Equation 
(133) for i equal to NS into Equation (101) with k subscript 
added to the blade root state vector yields 



Application of zero torque condition at pitch bearing: The 

condition of zero torque at the blade pitch bearing is achieved 
if the Laplace shifted (subscripted) form of Equation (104) which 
is 



0 


is satisfied. Insertion of Equation (133) with i replaced by N^ ea 
into the above expression and noting that the fourth element in 

each h 

l k,n 

equation 


CN fea> 


m 


is zero in value the pitch bearing discontinuity 



Application of control torque condition for coupling with 
swashplate or gyroscope control systems: The equation represent- 

ing the coupling between the blade and the swashplate motion 
when aerodynamics are allowed can be obtained by first considering 
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the modifications necessary to Equations (107) , (108) , and (109) , 

Shifting of the Laplace transform variable of Equation (107) 
yields 


k 1 (AT, ) - a m Y * 

m km mm 


- (N c t > 


00 i£(<f> -<J> ) 

+ I (1 + d m S £ )Y m e m mP w £ (s-ikn) = 0 (138) 

where 

Y* = a m [l + t (s-ikfif| 
m m |_ m 1 

Equations (108) and (109) can be extended to the subscripted forms 



and 


- <N c t > 
<+xk> m 




m 






(N 


H, 


c. t 


m 


) 


in which the Kb 


) 

terms have been dropped since the 

f- 1 (N c t } 

third element in each array is zero in value and the j c^. n V 


c l (N ct. 

k,nj 


m 


terms have been neglected since the pitch bearing should occur at 
or inboard of the control torque application point. The substi- 
tution of these last two equations into Equation (138) yields 
the control torque discontinuity equation in the form 
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For a gyroscope control system the last term in Equation (139) 
is replaced with a summation from £=-l to +1 of the term 

*1 (s - iks) - 

Application of flapwise moment condition at flap hinge: The 

condition of zero flapwise moment at the flap hinge is achieved 
if the Laplace shifted (subscripted) form of Equation (114) 
which is 


h>) 


(N flap> 


is satisfied. Insertion of Equation (133) with i replaced by 
Nfiap into the above expression and noting that the eleventh 


element in each •IcT, 


( k ,njm 

flapwise discontinuity equation 


(N fl ) 

v is zero in value results in the 


h'LL 


n=-N, 


^ N f lap^ 


(Tip) 


k , nl m 


(N flap> 


*c . t . 



In conjunction with the swashplate equations of motion. 
Equations (136), (137), (139), and (140) provide the necessary 

relations to obtain the solution when aerodynamic coupling is 
included. The state vector at any section can be obtained by use 
of Equation (133) after the values for tip Unknowns have been 
obtained. 
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Development of the Final Matrix Equations Governing 
Dynamic Forced Response Behavior of Rotor Blades 
Coupled Through a Swashplate 

The basic governing equations have been developed which are 
necessary to represent the dynamic forced response behavior of 
a rotor blade system coupled with a swashplate or gyroscope 
control system. These equations, however, have not been com- 
bined to construct the final form of the matrix Equations which 
must be solved in order to obtain blade tip and control system 
unknowns. To accomplish this, the swashplate governing equa- 
tions will be reduced to a coefficient form suitable for inclu- 
sion with an array coefficient form of the blade representation 
into a final matrix form. In addition, a similar coefficient 
form of the gyroscope equations will be constructed to replace 
the swashplate representation if this type of control is desired. 

General coefficient representation of swashplate governing 
equations of motion .- The swashplate governing equation of 
motion, Equation (71) , can be written in the form 


m(s 2 - 2i*,fis - £ 2 $2 2 ) + 
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(141) 


m=l 


m 


where 6^ is the Kronecker delta function defined such that 


6^—1 for i = j 
3 


and 


6^ = 0 for i ^ j 


Also, in obtaining Equation (141) , the force in the control rods, 
P m (s) , has been replaced by the control torques, (AT)^, divided 

by the respective control rod attachment point offset. 


On observation of Equation (141) , it appears that there would 
be more unknowns than equations because of the shifting of the 
Laplace transform variable. There is, however, an implicit 
relationship between w (s) and w (s-inft) . In this analysis, more 

X/ Xj 

equations are constructed from the basic ones by additional shift- 
ing of the Laplace transform variable such that s is replaced by 
(s-ikft) throughout Equation (141) . The end result in general form 
is 
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( 142 ) 




+ I (1 - b.S,)2 KC „ . + 


il+n 2 

R sp 


L + S 2 KC. 

6 JH-n ^1+n 



Note that (AT, ) which appears in Equation (142) is the same term 
k in 

which appears in Equation (138) and defined following that equa- 
tion. 
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Final matrix form of equation for determination of unknown 
quantities . - The complete set" of equations representing the 
coupled rotor-swashplate system can be written for the k denoted 
harmonic relative to that specified by the value used for the 
Laplace transform variable in the matrix form 


where 


n L Nf [Vn]{5 n }* + K}‘ " W 

[fk,n] 


k = 0, ±1, . . . ±N , 


(145) 


is representative of the terms multiplying the 


unknown quantities 


{%} 


and 


consists of the forcing function 


terms for all blades and swashplate. All the terms involved may be 
obtained from Equations (136) , (137) , (139) , (140) and (142) . 


The 


^k ; n 


array for a three-bladed rotor can be represented, 


utilizing integers for blade subscripts, as 
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8. 
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(146) 


where, if the number of spacial harmonics retained in the swash' 

plate is limited by letting the summation over £ range from 

-N to N , for N equal 1 the swashplate impedances from 
max max ' max ^ r 

Equation (142) are represented by 




(147) 
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and swashplate coupling to the blade, m, is given by 
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The blade root condition and discontinuity equations provide the 
generalized impedances of the blade, m, as 
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It should be noted that n j^* *"* / appearing in Equations (148) 

and (149) has a non-zero element only when k is equal to n. 

The column of unknowns, | q n j , represents the blade and 
swashplate unknowns in the form of 




n 
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where in general. 




w_ N (s-infi) 
"max 


w_ (s-inft) 
wq ( s-inn) 
wi (s-inH) 


w N (s-inft) 
max 


} and 


{P n }m = 



(152) 


The forcing function in Equation (145) can be expressed in 
the form 
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(153) 



\ j 


where 



* 

consists of 


(2N +1) 

max 


zero elements and 



(154) 


Modifications of final matrix representation for gyroscope 
control" system. - If a gyro system is used, the Laplace transform 
variable s can be shifted by ikft in Equations (20), (21) and (22) 
and manipulated similar to the swashplate governing equation, 
Equation (71) , such that the gyro governing equations can be 
written in the form similar to Equation (142) as 
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where 


w (s-inft) = (f> (s-inft) 
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K b k = CK’ + bfiC'i + C' (s-ikfi) ] . 

From Equation (155) , for a gyro instead of a swashplate, 
Equation (147) would be replaced with 


[f k ' n ] 



(156) 


(157) 


(158) 


87 



where use has been made of the fact that XG^'^ and VG°are zero. 
In addition, Equation (148) would be replaced with 




(159) 


The control torque discontinuity equation for the gyro representa- 
tion, Equation (113) , would result in the replacement of Equation 



(160) 


Having reduced the analysis to the matrix form of Equation 
(145) , the method of solution and wake iteration scheme can now 
be presented. 


Solution Scheme For Obtaining Forced Response 
Unknowns and Wake Iteration Procedure 

In the previous analytical section, the necessary expressions 
have been obtained to enable solution for the dynamic response 
variables for a helicopter rotor system coupled with a swashplate 
or gyroscope control system. This section is concerned with the 
basic solution to obtain values for the swashplate and blade tip 
unknowns, harmonically, and subsequent determination of the state 
vector variables and with the extension of this basic solution 
scheme to incorporate the effects of a deformed free wake on the 
dynamic response variable involved. 
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Solution scheme to obtain swashplate and blade tip unknowns . - 
In the development of analysis , use has been made of Laplace 
transform variables such as (s-inft) where the integer, n, repre- 
sents the harmonic shift relative to the particular harmonic of 
interest which has previously been denoted as Np. To solve for 
the Np/rev. swashplate and blade tip unknowns, the value of s in 
the analytical expressions necessary to obtain the terms involved 
in Equation (145) is taken as iNpft. Thus, these terms will be 
dependent on the value of Np in addition to the values of k and/ 

or n. Consideration of Equation (145) for -N^-k-N^ results in a 

matrix equation of the form exemplified by taking N^=l which is 
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and when multiplied by the inverse of the square matrix yields the 
solution for the unknowns as 
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arrays consist of Laplace transform operators and 

time dependent constants and the column arrays represent harmonics 
above, below, and at the harmonic of interest depending on the 
value of the subscript as referred to in the section pertaining 
to interpretation of the use of shifted Laplace transform vari- 
ables. Thus, for a given value of Np, N f , and N max , Equation (161) 

can be constructed and solved by Equation (162) for the (Np+l)/rev. 


where the 


T, 

_ k,n_ 



Np/rev., and (Np-1) /rev. blade tip and swashplate unknowns. 


Determination of the total state vector variable distribu- 
tions? - Once the unknowns are obtained, by Equation (162) for s 
equal to iNpft, Equation (133) is used with k equal to zero to 
determine the Np/rev. state vector variable harmonic coeffi- 
cients at the inboard end of each blade section. For a particular 

state vector variable, for example ,Np, which is the Np/rev. 

z 

harmonic coefficient of the shear force in the local z direction 
inboard of the ith section, is obtained in the form 

,Np = Ajg 1 + i . (163) 


As previously noted the state vector variables can be ex- 
pressed in a Fourier series form such that for example 
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Taking advantage of the fact that -n/rev. harmonic coefficients 
are complex conjugates of the +n/rev. harmonic coefficients, sub- 
stitution of Equation (163) into the last expression yields 
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where B 0 =0. 

This can also be written in the form 
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where 
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and 
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arctan 


(167) 
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If represents the azimuthal position of the blade relative 
to its position at t = 0, Equation (165) can be written as 
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(168) 


which represents the real value of the shear force acting perpen- 
dicular to the chord at ^ for section (i) . 


Other time dependent state vector variables can be obtained, 
utilizing the same representation as in Equation (168) where the 
harmonic coefficients used would correspond to the variable of 
interest. To obtain the total state vector variables as a func- 
tion of blade radius and azimuthal position the harmonic coeffi- 
cients for' each state vector variable must be determined for 
successively increasing values of N^ (beginning with zero) until 

the magnitude of the harmonic coefficients can be considered to 
be negligible compared to the lower harmonic coefficients. The 
harmonic coefficients for each state vector variable can then 
be utilized in the form of Equation (168) with the upper limit 
on the summation being determined by the highest non-negligible 
harmonic to obtain the radial and azimuthal distribution of the 
state variable of interest. 

Solution scheme for inclusion of deformed free-wake effects 
on dynamic response.- The inclusion of the effects of a free 
wake utilizes the previously discussed solution scheme used 
when free-wake effects are not included as an integral part of 
an iterative procedure which allows the blade response and free- 
wake to interact. 


The free-wake effect is initially taken into account by using 
the wake-induced velocity influence coefficient matrix (a-matrix) 
and the bound circulation matrix (r-matrix) generated by the free- 
wake program to compute the wake induced velocities at the blade 
aerodynamic section azimuthal locations. With the initial induced 
velocity distribution known, these values are introduced into the 
analysis for proper determination of angles of attack and aero- 
dynamic loads and moments at azimuthal positions for each aero- 
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dynamic section. This, information allows the construction of the 
aerodynamic transfer matrices necessary for the determination of 
the blade dynamic response in terms of a summation of the steady 
and harmonic response. 

From the dynamic response , in particular shear force 
quantities, a radial and azimuthal lift distribution acting on the 
blades can be obtained. This distribution includes the effects 
due to blade response which in turn is a function of the initial 
velocity field. To facilitate the determination of lift values, 
the blade sections with aerodynamics representation are taken to 
consist of only aerodynamic characteristics such that a differ- 
ence in the shear forces acting on the inboard end of the aero- 
dynamic blade section and inboard end of the previous outboard 
non-aerodynamic blade section can be utilized to obtain the 
actual lift forces acting on the aerodynamic blade section. For 
consistency in shear force vector orientation to insure proper 
shear difference determination both the aerodynamic and the pre- 
vious outboard non-aerodynamic blade section must have the same 
local coordinate system (i.e., the values for $, 0, ¥ must be 
identical for both sections) . 

The actual shear forces applied to an aerodynamic blade 
section, i, at the azimuthal position, ik , due to aerodynamic and 

blade response effects on this section can be expressed as 


W 


(i) 
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where the Ci~l). superscript denotes the previous outboard non-aero- 
dynamic blade section. The shear force, V z ($ ^ which has been 

defined by Equation (168) is the force acting on the outboard end 
of the (i+1) th section in the local blade coordinate system z 
direction (perpendicular to the blade section chordline not per 

turbed by cyclic pitch) . The shear force, -V (4^) ^ , which is 

defined by a similar form of Equation (168) is the force acting on 
the outboard end of the (i+1) th section in the local blade co- 
ordinate system -y direction Tpa^a-llel to the blade section chord- 
line not perturbed by cyclic pitch) . It should be noted that the 
direction of the local blade coordinate system z axis does not 
coincide with the direction of lift acting on the blade. 
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The lift distribution acting on the i th aerodynamic blade 
section at the kth. azimuthal position with variable induced velo- 
city and blade response effects included can be defined as 



= A 


VV 
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cos a 
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sin 



( 170 ) 


where the angle, a, , is the angle between the section chord line 
(zero cyclic pitch; and the resultant velocity vector acting on 
the ith aerodynamic blade station at the k th azimuthal position. 
This angle is defined approximately by the expression 


(i) 


~ a i ,k a c 


k 


where a is the cyclic pitch, defined positive for airfoil lead- 
c k 

ing edge rotated upward, and a. . is the angle of attack of air- 

X ,K 

foil including cyclic pitch at section, i, and azimuthal position 
denoted by i|). determined during the previous blade dynamic re- 
sponse calculations. Thus, in order to determine the lift distri- 
bution as in Equation (170) it is necessary to consider several 
harmonic coefficients of the shear forces for use in the form of 
Equation (168) and subsequently in Equations (169) . 


A new nondimensional circulation distribution, , pre- 

viously termed r-matrix, is obtained by use of the expression 


(i) ~ T (A) / _ ( i) 

k ~ L k / p V k 
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c ± ) 

where 1 is the effective length of the ith aerodynamic section 
and is amplitude of the resultant velocity vector determined 


during the previous blade dynamic response calculations acting on 
the ith aerodynamic section at the k th azimuthal location. In the 
aerodynamic representation the aerodynamic loading is applied at 
a point but the loading values represent that due to an effective 

spanwise aerodynamic length, L ^ , over which length the aerody- 
namic characteristics such as angle of attack and velocities are 
identical. This method of determining the new nondimensional cir- 
culation distribution is an approximation since both the angle of 
attack and velocity distributions used are those obtained during 
the previous dynamic response calculations (iteration) since the 
exact values are unobtainable. 
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The new is combined with the old [a] -matrix from the 

free-wake program tc obtain a new nondimensional induced velocity 
distribution which when multiplied by QR results in a new induced 
velocity distribution reflecting the coupling of the wake to the 
blade motions. This distribution can then be reinserted into the 
aerodynamic analysis and a new lift, circulation, and induced 
velocity distribution obtained. This process can be repeated 
until the desired degree of accuracy (or repeatability) is estab- 
lished at which time additional higher harmonics, not required 
during the wake iteration procedure, may be calculated. This 
section concludes the discussion of the analytical representation 
required to investigate the effects of a nonuniform swashplate on 
the dynamic response, and control loads of a helicopter. 


APPLICATION OF COMPUTER PROGRAM 


The computer program developed on the basis of the analyti- 
cal study was utilized to obtain the blade dynamic response of a 
four bladed H-34 rotor system model with various collective and 
cyclic control system stiffnesses. The results obtained showed 
only minor effects on blade dynamic response due to changes in 
the support system stiffness characteristics since the H-34 
system is stiff and restricts interblade coupling which would 
alter the natural frequencies and result in additional frequen- 
cies as was seen in the investigation of the much softer 0H-6A 
support system investigated in reference 2. The use of artifi- 
cial elastic axis and mass offsets originally used to enhance 
coupling, instead of the zero values the H-34 blade normally 
has, resulted in steady and harmonic pitch horn loads altered 
significantly from those obtained from experimental results. 

This section will be concerned with the flight condition and 
swashplate control system configurations for which the computer 
program was executed and some of the results obtained. It 
should be noted that the SI system of units will be used for 
purposes of discussion although the British system of units were 
employed in the actual program execution. 

The fundamental H-34 rotor system information was obtained 
from reference 3. The model of the blade (lumped parameter model) 
was based on the blade model employed for reference 1 with some 
modification to reduce the number of stations representing the 
mass, elastic, and geometric blade characteristics. This data 
included elastic axis and mass offsets which were originally 
thought to have little, if any, significant effect on the results. 
The final blade data employed is presented in Appendix A in tabu- 
lar form. The aerodynamic effects on the blade were included by 
using eight aerodynamic stations with effective aerodynamic 
lengths to represent the NACA 0012 airfoil which begins at 1.397 
meters from the axis of rotation and ends at the blade tip. The 
effective aerodynamic length of the outboard aerodynamic section 
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was modified to take into account tip losses. The blades had a 
uniform blade twist of -.139626 radians (outboard sections rotated 
leading edge down compared to inboard sections) over the blade 
airfoil length. On the basis of measurements reported in 
reference 3 the blade steady coning angle used was -.07854 radians 
(blade coned upward) and the lag angle used was -.10472 radians 
(blade lagged aft). The pitch bearing (feathering), lead-lag 
hinge (sweep) , and flap hinge were located .3048 meters from the 
rotor hub centerline and oriented going inboard to outboard in 
the order; flap, lead-lag, and pitch. 

The control rods were attached to the blade shaft by the 
pitch horn ahead of the blade pitch axis a distance of .20318 
meters. Attempts to. acquire information concerning the breakdown 
of the H-34 swashplate support system stiffnesses were not 
successful. Reference 3 provided a value of approximately 
78,638 Newton-meters per radian for the control system stiffness. 
Using this value as a basis, a collective system stiffness of 
32,091 Newton-meters per radian and lateral and longitudinal 
cyclic control system stiffnesses of 63,307 Newton-meters per 
radian were chosen as the reference uniform stiffness character- 
istics of the control system. Torsional response calculations 
were carried out for five flight conditions consisting of one 
cruise speed flight, y=.172; one high speed flight, y=. 291 ; one 
right turn, y=.235 and 1.34g loading? one left turn, y=.244 and 
1.35g loading; end for one pullup, u=. 223 and a 1.4g loading. 

The values for the advance ratio, y, given above are based only 
on the ratio of forward speed to rotational tip speed and do not 
include induced effects. Of these five flight conditions the 
first (cruise flight) was not based directly on a flight case in 
reference 3 since the desired advance ratio was not available 
whereas the remaining four are based in the order given on the 
experimental flights referred to as Flight 16, Flight 39, Flight 
40, and Flight 89 (Data Table 123), respectively, in reference 3. 
From these experimental conditions and ones similar to the cruise 
flight condition, values for the forward speed, rotational speed, 
air density, shaft tilt angles, collective pitch angle, and cyclic 
pitch angles were obtained. 

The program was initially executed for the five flight con- 
ditions with the reference uniform control system stiffness and 
with uniform downwash velocities based on Gessow and Myers (see 
reference 4) to ascertain if the necessary thrust was being 
obtained and also to obtain the steady and one per rev flap motions 
necessary for the wake geometry program developed in reference 1. 
The resultant steady thrust values per blade obtained from these 
runs were much less than the weight of the aircraft. To increase 
the thrust for each case, the values of lateral and longitudinal 
cyclic pitch angles were fixed slightly decreased from the orig- 
inal values while maintaining the same azimuthal position of 
maximum negative cyclic (defined as 6 radians ahead of the posi- 
tion ip=0 degrees in our convention or ip-90 degrees in the standard 
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convention) and the collective pitch increased. For each flight 
the effects of an increase in collective pitch angle increased 
the thrust in essentially a linear fashion (straight line curve) 
such that the necessary collective pitch values for the desired 
thrust could be obtained by interpolation. The criterion for 
the proper thrust per blade in forward flight was that the value 
for thrust should be between 12,788.6 Newton's (average weight 
of helicopter divided by the number of blades, four for the H-34 
helicopter) and 12,010.2 Newton's (the first value minus the 
blade weight) since the program internally considers the blade 
weight effects. For the maneuver type conditions the g factor 
was applied to the upper and lower limits to obtain the desired 
thrust range. The resulting collective and cyclic pitch input 
values obtained for each flight condition to yield the required 
thrust and, in addition, the uniform downwash velocities and the 
resultant thrust per blade are presented in Table I. The method 
utilized to obtain cyclic and collective pitch values represents 
only a trimming of control variables in regards to thrust values 
and does not include criteria related to roll and pitch moment 
values . 


The sets of control system stiffness values utilized by the 
program with various flight conditions were chosen to allow a 
rough determination of the general effect of changes in collec- 
tive stiffness, cyclic stiffness, or the control rod stiffness 
on the blade dynamic response behavior in addition to veri- 
fication of correct program operation. The parameters which 
were previously presented, necessary to define the lateral 
cyclic, longitudinal cyclic and collective control system stiff- 
ness, are the linear spring stiffnesses; , K and k^; the 

linear offsets; b., d / and a ; the location angles; <}> and x-; 

3 m m m 3 

and the integer controls; N, and N . With the assumption 

jO 6 • S • 

of four blades equally spaced and the offsets b^ and d^ equal to 

zero, the collective stiffness, lateral cyclic stiffness, and 
longitudinal cyclic stiffness can be expressed in the form*: 


K _ . = a 2 / 

coll. m ' 


N 

0 S 

^ + K b / }' k : + N b / K 

j = l J 


K 


lat. 



J k .cos 2 x • 
j=l ^ *3 


/ 


K, 

long. 



N 


e. s . 


I 

j=l 


k .sin 2 x • 
3 3 


/ 



In 1 

e.s. 


2 + 

I ^ 1 cos 2 x 1 

j=l 3 3 J 

/k 

m 


In 1 

e.s. 


2 + 

J k.sin 2 x- 

/k 
' m 


*The derivation of these equations is available in reference T 
which is the documentation pertaining to this computer program. 
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These expressions, assuming three equally spaced elastic supports 

2 TT 4 IT 

(i.e. N q s = 3, X-l - 0# X 2 ~ — » and x 3 = y), were used to 

obtain the control system collective, lateral cyclic, and longi- 
tudinal cyclic stiffnesses corresponding to given sets of input 
values for the variables associated with the linear spring-damper 
units involved in the representation of the swashplate control 
system. The input spring-damper stiffness values which were 
used and the resulting collective, lateral cyclic, and longi- 
tudinal cyclic control system stiffnesses are presented in 
Table II. 

Wake iteration program runs were executed for various 

combinations of flight conditions and control system stiffness 

configurations using the steady, one per rev., and two per rev. 

dynamic responses to determine the induced velocity distribution 

for each subsequent iteration. All of these runs are listed by 

flight conditions in Table III with information as to the rotor 

rotational speed, (radians/second) ; forward flight speed, 

V (meters/second) ; air density, p (kilograms/cubic meter) ; and 

shaft tilt angle, a (radians above the position of the shaft 

s 

parallel to the direction of forward flight) . The run code con- 
sists of a numeric character signifying the flight condition 
involved and an alphabetic character denoting the control system 
stiffness configuration used. The resulting steady blade root; 
pitch angle (positive leading edge rotated upward) , coning angle 
(positive blade tip rotated downward) , lead-lag angle (positive 
blade tip rotated forward) , and thrust per blade are also 
included. 

The blade root angles are obtained by adding the resulting 

blade root perturbation angles; $ , <j> , and <j> ; in respective 

x y z 

order to the collective pitch, coning, and lead-lag input values 
corresponding to the flight condition involved. The steady 
blade root pitch, coning, and lead-lag values in Table III are 
their respective values just outboard of the feathering bearing, 
flap hinge, and lead-lag hinge, respectively. The changes in 
the steady blade root coning and lead- lag angles from their 
respective input values are more pronounced than the changes in 
the steady blade root pitch angles from their input values. 

This is due to the fact that the input coning and lead-lag angles 
are only approximate and define the orientation of the blade 
spanwise axis in the local coordinate system (reference system 
for determination of response variables) as if the blade was 
rigidly rotated about the flap and lead-lag hinge axes. When 
the blade is allowed to respond to the applied forces and 
moments, as in a dynamic response calculation, the blade may 
deflect significantly from the reference positions at each 
radial section including blade root position. The change in 
steady root pitch is directly related to the control system 
collective stiffness which in the runs executed was sufficiently 
high such that only a very small steady pitch angle perturbation 
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TABLE III 

PROGRAM RUNS EXECUTED 


Run 

Flight 

Variable 

Steady Root Values 

Pitch 

(rad) 

Coning 

(rad) 

Lead- lag 
(rad) 

Thrust 

(N) 

1-C 

fi=22 . 4100 

V=55 . 5993 

p=l. 09775 

a =1.4661 
s 

0.29943 

-0.08112 

-0.13617 

12,082.7 

1-D 

0.29919 

-0.08126 

-0.13615 

12,108.1 

1-1 

0.30203 

-0.08369 

-0.13840 

12,428.3 

1-J 

0.29910 

-0.08151 

-0.13572 

12,147.6 

to 

1 

a 

0=22.4000 
V=32. 9581 
P=l. 08229 
a s =1.5208 

0.25773 

-0.08147 

-0.10053 

12,112.5 

2-D 

0.25754 

-0.08145 

-0.10003 

12,110.7 

2-E 

0.25672 

-0.08070 

-0.09983 

12,012.4 

2-F 

0.25726 

-0.08119 

-0.09997 

12,077.8 

2-G 

0.25495 

-0.07952 

-0.09814 

11,862.5 

2-H 

0.25548 

-0.08001 

-0.09828 

11,927.0 

H 

CNJ 

0.26019 

-0.08367 

-0.10123 

12,401.6 

2-J 

0.25748 

-0.08159 

-0.09943 

12,132.5 

3-A 

0=22. 8290 
V=45. 8179 
p=l. 08229 
a s =1.5202 

0.29767 

-0.10374 

-0.11018 

15,795.2 

3-B 

0.29868 

-0.10381 

-0.11170 

15,790.3 

4-A 

0=22.5150 
V=46. 8479 
p=l. 08229 
a =1. 5179 

0.29767 

-0. 10463 

-0.10809 

15,782.7 

4-B 

i 

0.29868 

-0.10467 

-0.11031 

15,773.8 

5-A 

0=23.0380 

V=43. 7602 

p=l . 08229 

a =1.6545 
s 

0.25697 

! 

-0.09688 

-0.15158 

15,138.2 
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was necessary to counteract the blade torque at the control 
torque application point in order to have zero torque at the 
pitch bearing. 

Comparison of blade root thrust values presented in Table I 
(thrust trim runs) and Table III (wake iteration runs) provides 
information as to the effect of inclusion of free-wake effects 
in the induced velocity distribution determination on the blade 
root thrust value. in the high speed forward flight condition, 
the thrust value dropped from 12,699.7 Newtons with a uniform 
induced velocity field to 12,082.7 Newtons with the free-wake 
generated induced velocity field (^5% drop) where the same 
control system stiffness, stiffness case C, was employed. In 
the cruise speed forward flight condition, the thrust value 
dropped from 12,481.7 Newtons to 12,112.5 Newtons (^3% drop). 

In a similar manner, using thrust values from Table III corre- 
sponding to the control system stiffness case B which has 
essentially the same collective control system stiffness as 
that in control system stiffness case C, the effect of including 
the free-wake effects in the maneuver flight conditions resulted 
in; a drop in the thrust from 16,792.0 Newtons to 15,790.3 
Newtons (^6% drop) for the right turn maneuver, a drop in the 
thrust from 16,769.8 Newtons to 15,773.8 Newtons (^6% drop) for 
the left turn maneuver, and a drop in the thrust from 17,543.8 
Newtons to 15,138.2 Newtons (^14% drop) for the pullup maneuver. 
Thus, the inclusion of free-wake effects in the blade response 
calculations significantly reduced the thrust per blade compared 
to that obtained with blade response calculations utilizing a 
uniform downwash distribution. However, the above percentages 
are dependent upon the input uniform induced velocity value 
used in each flight condition to determine the base thrust 
value. 

The independent effects of the control system; collective 
stiffness, lateral cyclic stiffness, and longitudinal cyclic 
stiffness; on the behavior of the four steady root variables 
presented in Table III cannot be ascertained in detail since 
the number of control system stiffness configurations for any 
one flight case is not sufficient to separate the cause and 
effect relationships involved. However, some general behavior 
characteristics of the steady root variables can be noted. 

Except for a few exceptions, the absolute values of the four 
steady root variables behave in the same manner as a result of 
a change in control system stiffness configuration. That is, 
for a given change in control system stiffness parameters all 
of the variables either increased or decreased in absolute 
value. Control system collective stiffness changes appear to 
be the dominant contributing factor to changes in the values 
for the steady root variables as would be expected. This 
conclusion is borne out to some extent by noting that of the 
runs pertaining to the cruise forward flight condition the 
highest absolute values for the steady root variables were 
obtained with the stiffest stiffness configuration, stiffness 
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case I, where the control system collective stiffness was 
60,935.7 Newtons-meters/radian and the lowest absolute values 
for the steady root variables were obtained for stiffness 
case G in which the control system collective stiffness was 
27,669.9 Newton-meters/radian. 

The results of the program runs were scanned to verify that 
the analysis operational in the computer program was performing 
correctly. The dynamic response variables; blade slopes, blade 
displacements, blade shears, blade moments, pitch horn loads 
(control loads), and swashplate control system deflections; 
resulting from each run in harmonic form were observed to 
satisfy all boundary conditions and discontinuity conditions 
employed by the program. Due to the magnitude of the values 
used in the control system stiffness configurations only minor 
changes resulted in the dynamic response variables as a result 
of changes in control system stiffness characteristics. Some 
of the harmonic results were employed to generate azimuthal 
distributions for comparison with experimental results obtained 
in reference 3. With the exception of pitch horn load distri- 
butions and torsional moment (blade torque) distributions the 
results compared favorably with the experimental results 
considering that the values for collective and cyclic pitch 
were different than those measured in reference 3 and that only 
the steady, one per rev. and two per rev', shear responses 
were utilized to obtain the free-wake generated variable in- 
duced velocity fields. 

The large discrepancy between the theoretically predicted 
and experimentally obtained pitch horn loading azimuthal distri- 
butions was investigated to ascertain its probable cause. This 
discrepancy in predicted and experimental pitch horn loadings 
for the high speed forward flight condition is depicted in 
Figure 10 , which represents the pitch horn loading azimuthal 
distribution obtained from program run 1-C and that obtained 
from reference 3. The sign convention for the pitch horn load- 
ing was based on assuming the pitch horn load to be positive 
when its effect is counteracting a positive blade torque (posi- 
tive for moment tending to rotate the blade leading edge upward) . 
Thus, for the H-34 helicopter rotor system in which the control 
rod is attached to the pitch horn ahead of the pitch axis, a 
positive pitch horn load represents a downward force applied to 
the pitch horn with the control rod being in tension. The 
azimuthal variable utilized in Figure 10 corresponds to the 
standard definition of azimuthal location (i.e. ^g^andard 

zero when blade is directly aft) and not the azimuthal reference 
variable utilized in the analysis. This required the theoretical 
pitch horn loading to be determined relative to the reference 
variable, ip, used in the analysis and then plotted at ^standard 

equal to ip + tt/2. The values for the pitch horn loading har- 
monics were obtained by dividing the (AT) variable outputted 
in the solution vector for each harmonic by the distance, a^, 
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Figure 10. Azimuthal pitch horn loading for 
the high speed flight 
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Figure 11. Harmonic content of pitch horn 
loading for high speed flight 



Figure 12. Phasing associated with the 

harmonic content of pitch horn 
loading for high speed flight 


where, depending on the blade model, the (AT) harmonic is equal 
or approximately equal to the blade torque harmonic of the 
section outboard of the section inboard of which the control 
torque is applied. 

On comparison of the theoretically and experimentally deter- 
mined pitch horn loading distributions it is obvious that a 
significant difference in the steady component occurs. The 
experimental results obtained from reference 3 were harmonically 
analyzed to obtain the harmonic amplitudes and associated phase 
angles for comparison with the corresponding theoretical values. 
The pitch horn loading harmonic content, amplitude and phasing, 
of the theoretical and experimental distributions shown in 
Figure. 10 is compared in Figures 11 and 12, The harmonic content 
shown in these figures has been determined relative to the 
azimuthal position defined by the reference variable, if), equal 
to zero. By taking the theoretical steady pitch horn to be that 
of the experimental data and using the theoretical harmonics, 
a third distribution (shifted) was obtained and added to 
Figure 10. The shifted theoretical curve is in better agreement 
with the experimental results. It can be concluded that a 
primary difference in the experimental and theoretical results 
is the magnitude and sign of the steady pitch loads. 


In an investigation to determine the cause of the large 
steady negative pitch horn loads, it was found that the mass 
and elastic offsets reported for the H-34 blade in reference 1 
and used in the present analysis were the significant contri- 
butors to large steady feathering torques which are directly 
related to steady pitch horn loadings. In the high speed 
flight condition the center of mass offsets of 0.00417 meters 
ahead of the neutral axis combined with perturbation quantities 
and centrifugal force components for all mass stations occurring 
outboard of the airfoil cutout was found to contribute approxi- 
mately 85 percent of the negative feathering torque acting at 
the control torque discontinuity application’ point. Roughly, 
an additional 10 percent of the negative feathering torque was 
due to the aerodynamic forces acting at quarter chord which 
was 0.002083 meters aft of the elastic axis. In addition, 
these artificial offsets may result in differences in the 
harmonic content of feathering torque and thereby pitch horn 
loadings both due to their involvement in the analytical rotor 
system representation for all dynamic response harmonics and 
due to interharmonic coupling in which the harmonics above and 
below the harmonic of interest can effect the results obtained 
for the harmonic of interest. Since these small offsets were 
not present in the physical system, it was concluded that these 
offsets were a major cause of the difference between the experi- 
mental and predicted results. This effect of the small offsets 
suggest that the dynamic response of a helicopter rotor system 
may be much more sensitive to the spanwise distribution of the 
center of mass and elastic axis locations than previously 
suspected. 
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Representative azimuthal plots of the predicted flapwise 
bending moment, chordwise bending moment, and torsional moment 
for the same high speed flight condition are compared with the 
measured results from reference 3 in Figures 13 thru 15. Since 
the H-34 blade system is highly uncoupled it is not expected 
that the neutral axis and mass offsets would provide a very 
significant effect of the high negative torsional loading on 
the bending moments . This is borne out in that the behavior of 
the predicted results seems to follow the azimuthal trend of 
the experimental data. The comparison of experimental and 
theoretical results would be further improved if additional 
harmonics were included in the interharmonic coupling repre- 
sentation (one/rev. above and below the harmonic of interest 
was used in program runs) and ir the free-wake induced velocity 
determination and if the artificial offsets were removed. 


CONCLUDING REMARKS AND RECOMMENDATIONS 


The program which has been developed for this investigation 
is a highly sophisticated theoretical representation of heli- 
copter rotor systems. The coded analysis includes the effect 
of an anisotropically mounted flexible swashplate or a gyroscope 
control system and the effects of a deformed free-wake. The 
analysis has the capability of predicting the response charac- 
teristics of helicopter rotor systems in both forward flight 
and steady-state maneuver conditions. The quasi-steady blade 
aerodynamics representation can utilize an induced velocity 
field which is either uniform, variably defined, or one which 
is generated on successive iterations based on information 
from a free-wake analysis. The programmed analysis includes 
aerodynamic interharmonic coupling and also structural inter- 
harmonic coupling thru the swashplate or gyroscope without 
the necessity of information concerning the rotor system 
natural frequencies or related mode shapes. The analytical 
method employed allows for the straightforward extension of 
the program to include unsteady phenomena such as transient 
maneuvers and gust loading effects. Also the program is in a 
form such that it could be modified to investigate the stability 
and vibration characteristics of helicopter rotor systems in 
forward flight or maneuvers. 

The blade dynamic response calculations which were carried 
out indicate that even with the high control system stiffnesses 
involved that the program provides an effective method of inves- 
tigating the dynamic response behavior of present or future 
helicopter rotor systems. Considering the extensive information 
that can be obtained by use of this program, the program is a 
very efficient tool for use in investigating the effects of 
various parameters on the dynamic characteristics of a heli- 
copter rotor system. A significant result of the execution of 
the program for the H-34 rotor system was the observation that 
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Figure 15. Torsional moments for the high 
speed flight 



there is a significant effect of small offsets of the center of 
mass from the blade neutral axis and of the blade neutral 
axis from the blade aerodynamic center on the pitch horn loads 
and torsional moments. This result suggests that the blade 
pitch horn control loads may be more sensitive to the spanwise 
distribution of center of mass, neutral axis, and aerodynamic 
center locations than previously thought. 

It is recommended that the program be extended to repre- 
sent the unsteady phenomena a helicopter system might encounter 
as well as altering the program to provide the option of calcu- 
lating the vibrational or stability characteristics of helicopter 
rotor systems. In addition, it is recommended that in order to 
more fully investigate the capabilities of the program and to 
show the effects of a nonuniform swashplate that the program be 
used to investigate the dynamic characteristics of a rotor 
system having a soft control system such as the 0H-6A. 
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APPENDIX B 


CONVERSION OF UNITS 

The following conversion relationships are those which were 
utilized to convert variable values originally in the British 
system of units to the corresponding values in the SI system of 
units for purposes of discussion in this report. 

1 foot = .3048006 meters 
1 pound = 4.448222 Newtons 
1 slug = 14.59388 kilograms 
1 foot/second = .3048006 meters/second 
1 foot-pound = 1.355821 Newton-meters 

1 pound/foot = 14.59388 Newtons/meter 

2 2 
1 pound-foot = .4132550 Newton-meters 

2 2 

1 slug-foot = 1.355821 kilogram-meters 

3 3 

1 slug/foot = 515.3750 kilograms /meter 

1 foot-pound/radian = 1.355821 Newton-meters/radian 

1 radian/ (foot-pound) = .7375607 radians/ (Newton-meter) 
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